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Preface 


The  objective  of  this  study  is  to  establish  the  capability  of  the 
Nielsen  store  separation  and  trajectory  program  to  predict  off-body 
store  loadings.  The  Nielsen  program  was  developed  in  the  early  1970's 
and  was  designed  specifically  to  handle  the  store  trajectory  problem. 

In  order  to  establish  the  program  capabilities,  a  computer  model 
of  an  advanced  fighter  design  was  generated.  The  specific  model  was 
chosen  because  it  had  been  extensively  studied  in  a  wind  tunnel  and 
aerodynamic  data  were  available  for  comparison. 

A  secondary  objective  is  to  consolidate,  as  much  as  possible,  the 
theory  required  to  understand  the  program  and  provide  the  aircraft 
computer  model  for  future  study.  Theory  notation  was  a  problem  because 
the  theory  was  developed  in  three  reports  with  the  same  symbols  used 
for  different  quantities  in  different  reports.  For  this  reason  a 
notation  section  is  included  for  those  paragraphs  where  clarification 
is  required. 

The  computer  program  is  not  included  in  this  report;  however,  the 
basic  program,  with  no  changes,  was  used  for  the  analysis. 

I  would  like  to  express  my  sincere  gratitude  for  the  assistance 
I  received  during  the  course  of  this  study.  My  primary  independent 
study  advisor.  Major  Michael  L.  Smith,  AFIT/ENY,  provided  the  advice 
and  guidance  which  allowed  me  to  complete  the  work.  I  owe  a  special 
debt  of  gratitude  to  Mr.  Richard  Dyer,  AFWAL/FIMM,  who  was  the  sponsor 
of  the  study.  I  also  owe  special  thanks  to  Mr.  Bart  Heath,  AFWAL/FIGC, 
who  provided  technical  guidance  on  the  use  of  the  computer  program. 
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b  fill  s 
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a  fuselage  and  store  angle  of  attack 

ct^  local  angle  of  attack  due  to  wing  twist  and  camber 
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e  polar  angle,  see  figure  2 

v,t  vectors  normal  to  and  tangential  to  the  crossflow  plane 

contour,  see  figure  2 

€,n,s  axial,  lateral,  and  vertical  displacements,  respectively, 

of  the  store  center  of  gravity  relative  to  the  carriage 
position  on  the  pylon  in  the  fuselage  coordinate  system; 
positive  forward,  to  the  right,  and  downward,  respectively 

P  angle  defined  by  equation  (26) 

<j>  complete  potential  for  a  body  defined  by  equation  (3) 

<l>e  potential  associated  with  equivalent  body  of  revolution 

. v  dihedral  angle  of  image  horseshoe  vortex 

<Pv  dihedral  angle  of  the  vth  vortex  on  the  exposed  wing  panel 
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crossflow  plane  potential  associated  with  the  contour  of 
a  general  body 


local  body  cross-sectional  area 

rate  of  change  of  body  cross-sectional  area 

2 

reference  area,  ird  /4 
time,  seconds 

incompressible  perturbation  velocities  in  x,y,z  directions 
of  figure  1;  compressible  perturbation  velocities  in  x,y,z 
directions  of  figure  1 

incompressible  externally  induced  perturbation  velocities 
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incompressible  perturbation  velocities  in  x,v,x  coordinate 
system 

incompressible  free-stream  velocity  components  in  x,v,x 
coordinate  system 

incompressible  radial  velocity  in  y-z  plane  of  figure  1 

incompressible  free-stream  velocity  components  in  y  and 
z  directions 

velocity  inside  engine  inlet 

compressible  velocities  in  store  body  coordinate  system, 
see  figure  10 

free-stream  velocity 

free-stream  velocity  vector  in  incompressible  space 

separated  store  free-stream  velocity 

coordinate  system  fixed  in  nose  of  axisymmetric  body  in 
incompressible  space,  see  figure  1;  coordinate  system 
fixed  in  nose  of  general  noncircular  body  in  incompres¬ 
sible  space,  see  figure  2;  coordinate  system  with  origin 
at  midspan  of  the  bound  leg  of  a  horseshoe  vortex  in 
incompressible  space 


angle  defined  by  equation  (38) 

sweep  angle  of  the  bound  leg  of  an  image  horseshoe  vortex 

sweep  angle  of  the  bound  leg  of  a  horseshoe  vortex  on  the 
exposed  wing  panel 

changes  in  yaw  and  pitch  angles,  respectively,  of  store  from 
time  =  0  values;  positive  nose  to  the  right  and  nose  up 
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An  advanced  fighter  design  is  modeled  using  the  Nielsen  store 
separation  and  trajectory  program  with  the  non-circular  fuselage  cross 
section  option.  The  theory  required  to  build  the  computer  model  is 
consolidated  for  better  understanding,  with  the  major  points  refer¬ 
enced  to  the  appropriate  report.  The  details  of  the  process  to  gener¬ 
ate  the  computer  model  are  described.  The  computer  model  of  the 
aircraft  and  ogive  store  is  used  to  predict  store  forces  and  moments 
which  are  compared  with  wind  tunnel  results  on  the  same  configuration. 
The  prediction  accuracy,  at  Mach  number  0.6,  at  a  store  position  of 
three  store  diameters  away  from  the  aircraft,  is  within  ten  percent 
of  experiment  in  the  region  of  interest,  which  is  the  area  under  that 
portion  of  the  fuselage  occupied  by  the  wing. 
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I.  INTRODUCTION 


The  analysis  of  what  happens  to  a  weapon  when  it  is  released  from 
an  aircraft  is  a  major  concern  in  weapons  development.  The  weapon  must 
separate  from  the  aircraft  without  structural  interference  and  it  must 
then  follow  a  predictable  path  to  impact.  The  prediction  of  the  sep¬ 
aration  and  the  initial  weapon  path  is  the  problem  to  be  addressed. 

The  solution  to  the  problem  requires  that  the  forces  on  the  weapon  be 
known.  These  forces,  and  the  equations  of  motion  predict  the  path  the 
weapon  will  follow. 

A  computer  prediction  method  was  developed  by  Nielsen  Engineering 
and  Research,  Inc.,  under  contract  to  the  United  States  Air  Force. 

The  work  was  done  during  the  period  1968  to  1972.  The  final  result  was 
a  method  for  predicting  the  six  degree-of-freedom  store  separation 
trajectory  at  speeds  up  to  the  critical  speed. 

The  procedure  was  developed  in  three  stages.  Stage  one  was  a 
three-degree-of-freedom  theory  (Ref  1).  Stage  two  was  a  six-degree-of- 
freedom  program  which  modified  the  three-degree-of-freedom  program 
(Ref  2).  This  second  stage  included  a  users'  manual  for  the  computer 
program  (Ref  3).  The  final  stage  extended  the  original  six-degree-of- 
freedom  approach  to  handle  more  general  configurations  (Ref  4,5). 

The  final  program  can  model  a  wing  with  twist  and  camber.  The 
wing  may  have  leading  and  trailing  edge  sweep,  and  dihedral.  These 
parameters  may  vary  along  the  span.  Pylons  and  racks  can  be  modeled 


although  the  connectors  between  store  and  rack  are  not  modeled.  Up  to 
ten  stores  may  be  included,  either  in  multiple  ejector  rack  (MER)  or 
triple  ejector  rack  (TER)  configuration.  The  store  is  modeled  as  a 
body  of  revolution  with  planar  or  cruciform  empennage.  The  store  is 
non-powered,  but  may  be  given  an  initial  velocity  and  orientation.  The 
fuselage  cross  section  can  be  noncircular,  but  it  must  have  a  vertical 
plane  of  symmetry.  Air  inlets  may  be  included  and  the  ratio  of  the 
velocity  in  the  inlet  to  the  free  stream  velocity  can  be  specified. 

The  aircraft  must  be  operating  at  constant  velocity  with  zero  yaw  angle 
and  constant  angle  of  attack.  The  computer  program  uses  a  constant, 
preselected,  atmospheric  density.  This  does  not  restrict  the  flight 
path  to  horizontal  only,  but  the  density  does  not  change. 

The  objective  of  this  report  is  to  validate  the  Nielsen  program 
by  applying  it  to  an  advanced  fighter  design.  The  store  to  be  separ¬ 
ated  is  a  generic  store  with  known  characteristics.  The  fighter  and 
store  have  been  aerodynamical ly  tested  in  a  wind  tunnel  and  these  data 
are  used  to  compare  to  the  computer  results. 

In  the  process  of  building  the  computer  model  it  was  noted  that 
the  theory  for  the  model  was  contained  in  three  different  reports. 

This  made  the  theory  very  difficult  to  follow.  As  a  consequence,  the 
theory  has  been  consolidated  with  each  section  referenced  to  the 
original  report.  Only  those  sections  essential  to  the  construction  of 
the  computer  model  are  included  in  detail. 
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The  details  of  the  aircraft,  store  and  wind  tunnel  test  conditions 
are  given  in  Appendix  A.  For  simplicity  and  data  matching,  only  the 
aircraft  and  two  stores  were  modeled.  No  pylons  or  racks  were  included. 
The  aircraft  is  an  advanced  fighter  design  with  a  non-circular  fuselage 
cross  section  including  flow-through  engine  inlets.  The  stores  were 
identical.  The  store  body  was  a  circular  cylinder  with  an  ogive  nose 
and  cruciform  empennage.  The  computer  models  were  used  to  generate 
force,  moment,  position,  and  relative  velocity  data  on  one  store  which 
was  placed  at  various  positions  under  the  fuselage.  The  trajectory 
option  was  not  used. 

At  Mach  numbers  0.6  and  0.9,  one  store  was  placed  approximately 
three  store  diameters  below  the  fuselage  centerline  and  moved  from  a 
position  under  the  fuselage  nose  to  the  tail.  At  Mach  number  0.6,  the 
store  was  placed  three  feet  to  the  right  of  centerline  and  moved  from 
nose  to  tail.  The  final  position  was  approximately  one  store  diameter 
below  the  centerline.  For  this  case  the  Mach  number  was  0.6  and  the 
position  was  varied  from  37.22  feet  behind  the  aircraft  nose  to  the 
tail.  One  store  remained  on  the  centerline  and  flush  with  the  fuse¬ 
lage. 

The  data  obtained  were  compared  with  wind  tunnel  values  on  the 
same  configuration.  The  primary  region  of  interest  for  store  position 
was  the  area  under  the  fuselage  in  the  range  occupied  by  the  wing.  At 
Mach  number  0.6  with  the  store  positioned  three  diameters  below  the 
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fuselage  centerline,  the  predicted  results  were  within  ten  percent  of 
experimental.  Accuracy  decreased  with  lateral  movement  under  the 
fuselage.  The  prediction  error  was  20  percent  or  greater  when  the 
store  was  moved  to  a  position  one  store  diameter  below  the  fuselage 
centerline. 


II.  THEORY-NIELSEN  PROGRAM 


The  three  principal  tasks  in  the  prediction  of  a  store  trajectory 
are:  first,  the  determination  of  the  nonuniform  flow  field  in  the 
neighborhood  of  the  ejected  store;  second,  the  determination  of  the 
forces  and  moments  on  the  store  in  this  flow  field;  and  third,  the 
integration  of  the  equations  of  motion  to  determine  the  store 
trajectory  (Ref  1:1). 

Compressibility  Correction  (Ref  1:4-5) 

Nomenclature: 

<j>  -  potential  function  in  compressible  space 

<PXX,  <Pyy,  *zz  -  32$/3X2  ,  32<j>/3y2  ,  32<|>/3Z2 

x,y,z  -  coordinate  system  in  compressible  space 
x',y‘ ,z*  -  coordinate  system  in  incompressible  space 

It  is  assumed  that  the  flow  field  for  compressible  subcritical 
flow  is  governed  by  a  perturbation  potential  $  which  satisfies 

U-H„!)+xx  +  +yy  +  *zz  ■  0  (!) 

The  flow  field  represents  a  wing-fuselage  combination  moving  at 
with  a  fuselage  angle  of  attack  of  og.  The  x,  y,  z  coordinate  system 
is  fixed  in  the  fuselage,  see  figure  1. 

The  wing  fuselage  combination  is  transformed  to  an  equivalent 
incompressible  one.  The  incompressible  flow  field  for  the  equivalent 


body  is  determined,  and  then  transformed  back  to  the  compressible  space 
The  following  transformation  is  used  from  the  compressible  space 


(x,y,z)  to  the  incompressible  space  (x',y',z'): 

x'  a  ~ — ~ »  y'=y  »  z'=z  (2) 

00 

Mathematical  Models  for  Determining  Flow  Field 

All  of  the  models  described  are  for  the  equivalent  incompressible 
configuration,  although  the  primed  notation  has  been  dropped.  This 
configuration  is  determined  by  the  program  from  the  input  of  the 
actual  aircraft  dimensions  (Ref  4:2). 

Fuselage  Potentials  (Ref  4:4)  The  model  of  bodies  with  non-cir¬ 
cular  cross  section  is  based  on  the  equivalence  rule  which  states  that 
the  following  conditions  hold  for  a  general  slender  body  (Ref  6:107-110) 

a.  Far  away  from  a  general  slender  body,  the  flow  becomes  axi- 
symmetric  and  equal  to  the  flow  around  the  equivalent  body  of  revolu¬ 
tion. 

b.  Near  the  general  slender  body,  the  flow  differs  from  that 
around  the  equivalent  body  of  revolution  by  a  two-dimensional  constant 
density  crossflow  part  that  makes  the  tangency  condition  at  the  actual 
body  surface  be  satisfied. 

The  outer  axisymmetric  flow  is  given  by  the  potential  4>e  which  is 
associated  with  an  equivalent  body  of  revolution.  This  equivalent  body 
of  revolution  is  axisymmetric,  with  the  same  cross  sectional  area  dis¬ 
tribution  along  its  longitudinal  axis  as  the  actual  body. 


The  inner  flow  is  represented  by  a  two-dimensional  potential  <|> 

2 

which  in  the  limit  (outer)  becomes  [Ux(x)/2TrVoo]Sl (x)ln  r.  The  total 
velocity  in  the  axial  direction  x  of  the  body  is  Ux(x).  The  rate  of 
change  of  the  cross  sectional  area  S  with  x  is  S'(x).  The  radial 
distance  from  the  body  centerline  is  r,  see  figure  2. 

A  solution  valid  for  the  entire  flow  field  is  given  by: 

U  (x) 

♦c(r,e)  =  <t>e(r)  +  *  (r,e)  -  S'(x)  In  r  (3) 

2  00 

The  potential  <|>e  and  $  are  determined  independently.  The 
potential  <|>e  of  the  equivalent  body  is  determined  from  flow  modeling 
of  a  body  of  revolution  at  zero  angle  of  attack  with  the  same  cross 
section  area  distribution  as  the  actual  body.  The  inner  potential  <{> 

2 

is  determined  in  the  crossflow  plane  and  satisfies  the  flow  tangency 
condition  on  the  actual  body  contour. 

Outer  Potential  -  4>c  (Ref  1:9-12) 

Nomenclature: 

aR  -  body  length 

Qk  -  strength  of  k  point  source 

V  V4*‘rv» 

r*-  r/*R 

U  -  axial  velocity 
U*  -  U/V 
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The  potential  4>e  itself  is  not  calculated,  although  the  pertur¬ 
bation  velocity  field  can  be  determined  by  knowing  the  strengths  and 
distribution  of  three-dimensional  point  sources  along  the  body  longi¬ 
tudinal  axis. 

The  stream  function  tp  and  the  axial  and  radial  velocity  components, 
U  and  Vr,  respectively,  are  given  by  (Ref  8): 


12  N  1  +  (x*-x*  ) 

=  _  r  n  * _  K 


<l<*(x*,r*)  *  fr*  -  l  Q  < 

k=l  K 


[(x*-xk*)2  +  r2]* 
U*(x*,r*)  =  A  78S  =  1  +  1  %* - — - — 


(4) 


k=rk  C(x*-xk*)2  +  r*z]3/2 


(5) 


V  *(x*  r*}  - _ L  ty *  r 

V  ,r  j  r*  dx*  i 


N  Qk*r* 

k=l  [(x*-xk*)2  +  r*2]  3/2 


(6) 
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Up  is  a  reference  length  which  will  be  taken  as  the  length  of  the  body 
(figure  1).  U  and  Vr  are  due  jointly  to  a  free  stream  velocity 
aligned  with  the  body  axis,  and  a  series  of  N  point  sources  distributed 
along  the  body  axis. 

The  quantities  Q^*  and  x^*  are  the  source  strength  and  location, 
respectively,  of  the  k**1  point  source  and  the  point  (x*,r*)  designates 
the  field  point  in  cylindrical  coordinates. 

The  flow  field  about  the  axi symmetric  body  is  obtained  by  imposing 
three  conditions  on  equations  (5)  and  (6).  These  velocity  equations 
are  used  to  calculate  the  source  strengths  rather  than  the  streamline 
equation  because  they  were  found  to  give  a  better  shape  representation. 
This  is  due  to  the  fact  that  they  require  both  body  ordinates  and 
surface  slopes  for  their  solution. 

The  first  condition  requires  the  flow  directions  at  (N-2)  points 
(x*,  r*),  j=l,...,(N-2),  be  specified  by  angles  e.  with  respect  to  the 

J  J  J 

positive  x*-axis: 

V  *(xT,z*) 

tan  b.  =  — — -J — 2—  for  j=l , . . . , (N-2 )  (7) 

u*(xj,zj) 

The  locations  (x*,r*)  where  the  flow  angles  B .  are  specified  lie  on 

J  J  J 

the  surface  of  the  axi symmetric  body  being  represented  and  the  values 
of  tan  8-  correspond  to  the  local  body  surface  slopes  dr*(x*,r*)/dx*. 

J  J  v 


Writing  equation  (7)  in  terms  of  equation  (5)  and  (6)  and  rearrang¬ 


ing  them  gives  a  set  of  linear  algebraic  equations  in  Qk*'s. 
N 


tan  e 


=  I  Q  * 
J  k=l  k 


r*  - ta"  ei(x*  -  V> 

2  3/2 

[(x*-xk*)  +  rj2] 


(8) 


The  second  condition  is  that  the  sum  of  all  the  source  strengths 
be  zero. 


N 

1  Qk*  =  o 

k=l  K 


(9) 


This  insures  that  the  surface  described  by  ^*(x*,r*)=0  will  be  a  closed 
surface. 

The  third  condition  is  the  existence  of  a  stagnation  point  at  the 
body  nose.  U*=0  at  x*=r*=0.  Substituting  this  into  equation  (5),  the 
third  condition  becomes 


(10) 


Since  the  sources  are  distributed  along  the  positive  x*-axis  only,  the 
forward  tip  of  the  body  is  positioned  at  the  origin  and  the  ip*=0  surface 
is  the  body  surface. 

When  the  N  source  positions,  xk*  are  selected,  equations  (8),  (9), 
and  (10)  comprise  a  set  of  N  linear  equations  in  N  unknown  Qk*.  Once 
Qk*'s  are  solved  for,  equations  (5)  and  (6)  can  be  used  to  determine  the 
perturbation  velocities  at  any  point  in  the  flow  field  around  the  body. 


Inner  Potential  -  <j>^  (Ref  4:4-18)  The  inner  potential  <j>  is 
composed  of  higher  order  singularities  given  by  polar  harmonics  and  a 
two-dimensional  source  term: 

MH  acosne  Uv(x) 

*  (r,e)  =  l  — - -  +  j-rj —  S'  (x)  In  r  (11) 

2  n=l  rn(e)  ^v- 

where  MH  is  the  number  of  polar  harmonics.  This  assumes  flow  symmetry 
about  the  vertical  x-z  plane.  Side  slip  is  not  considered.  Two- 

dimensional  polar  harmonics  are  described  in  greater  detail  (Ref  7: 
40-48,  55-59). 

Polar  harmonic  solutions  are  obtained  at  a  number  of  crossflow 
plane  stations  along  the  body  longitudinal,  x,  axis.  The  x  range  is 
the  range  over  which  the  local  flow  needs  to  be  calculated.  At  each 
crossflow  plane  station,  a  number  of  control  points  are  distributed  at 
equal  angular  spacings  in  0  on  the  body  contour  over  the  range,  0=0° 
to  0=180°,  see  figure  2. 

The  flow  tangency  condition  is  applied  at  the  control  points  re¬ 
sulting  in  a  finite  set  of  simultaneous  equations  in  terms  of  the 
unknown  polar  harmonic  coefficients. 

To  apply  the  flow  tangency  boundary  condition  a  coordinate  system 
is  set  up  at  each  control  point,  figure  2.  The  system  is  set  up  such 
that  v  is  in  the  direction  normal  to  the  contour  and  t  is  tangential 
to  the  contour.  Figure  2  shows  unit  vectors  associated  with  those 
directions.  It  is  assumed  (Ref  6)  that  the  shape  of  the  body  surface 
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can  be  expressed  as: 


F(x,v,t)  =  0  -£<x<0  (12) 

where  i  is  the  body  length.  The  unit  vector  normal  to  the  body  surface, 
not  the  contour,  at  the  point  is  given  as: 

®n  *  ^  <“> 

where 

5v+f  5x  (14> 

Let  the  free-stream  velocity,  which  does  not  have  to  be  uniform,  and 
the  perturbation  velocity  vector  be  denoted  ^.“^(x.v, t)  and  q^q^x, 
v , t ) ,  respectively. 

The  flow  tangency  boundary  condition  can  now  be  written  in  the 
(x , v , t )  system 

[Vjx.v.t)  +  qoo(x,v,T)]*in  =  0  (15) 

Expanding  V^,  qro  in  terms  of  their  components  in  the  x,v,  and  t 
directions  results  in 

^(x,v,t)  +  q  Jx.v.x)  =  (Ux  +  ux)ex  +  (Uv  +  uv)ey 

+  (U^  +  ut)eT  (16) 

Substituting  equations  (13),  (14),  and  (16)  into  the  boundary 
equation  (15)  gives 

(U  +  u)  —  +  (U  +  u  )#  +  (U  +  u  )  -—  =  0  (17) 

v  X  Xy  9X  '  v  v'  9v  v  T  i'  3t  ' 
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This  equation  represents  the  nonlinear  flow  tangency  condition  for  any 
body  situated  in  a  non-uniform  flow. 


Equation  (17)  will  now  be  simplified  with  no  loss  of  generality. 
The  unit  vectors  5n  and  5^  at  a  point  on  the  body  surface  lie  in  the 
same  longitudinal  plane,  see  figure  2.  Since  5^  is  perpendicular  by 
definition  to  e^,  then  in  is  also  at  a  right  angle  to  i^.  Therefore: 

V  eT  =  0  (18) 
Using  equations  (13)  and  (14),  equation  (18)  becomes: 


(—  5 
'  3X  x 


e 

t 


=  0 


(19) 


Since  F(x,v,x)  is  a  constant  and  equal  to  zero  on  the  body  surface,  it 
does  not  vary  with  respect  to  the  t  direction 

|£  =  0  (20) 


Taking  the  differential  of  F(x,v,x)=0 

dF(x,v,x)  =|^dx  +|^dv  +|^dx  =  0 

With  equation  (20)  this  becomes: 


=  .  if.  dv 
ax  av  3x 


(21) 


(22) 


Substituting  equations  (20)  and  (22)  into  (17)  gives  the  nonlinear  flow 
tangency  condition  that  must  be  satisfied  at  points  on  the  surface  of 
a  general  body. 


(23) 


jg  *C i 


[Ux(x)  +  ux(x,v,t)]  ^  =  Uv(x)  +  u  (x.v.t) 


The  problem  is  to  determine  ^  and  then  write  Uv  and  uy  in  terms 
of  y  and  z  components  along  the  e  and  r  directions. 

An  expression  for  dv/dx,  the  streamwise  body  slope  will  be  de¬ 
rived  in  accordance  with  the  principles  of  Ref  6,  see  figure  3.  The 
body  radius  r  and  its  derivatives  with  respect  to  e  must  be  single 
valued  and  continuous  for  a  given  value  of  e. 

Vector  v  is  perpendicular  to  and  t  tangential  to  the  body  contour 
at  x.  Let  Av  denote  the  change  measured  in  the  direction  of  v  of  the 
location  of  the  contour  when  going  in  a  longitudinal  plane  from  the 
cross  section  at  x  to  the  one  at  x-ax,  see  figure  2  for  definition  of 
x  direction.  The  slope  of  the  body  in  the  longitudinal  planr  is 


.  linl  ivM 

dx  AX 


This  body  slope  is  approximated  by 


dv(e)  _  Ar(e)cos 
dx  “  Ax 


4r<6>  =  rX-4X(9>  -  rx<e’ 


From  figure  3,  angle  <p  between  the  normal  v  and  the  radial  direction  r 
can  be  expressed  in  terms  of  the  contour  slope  m  and  the  polar  angle  e: 


<j»=m-0- 


The  body  contour  shape  m(0<m<27r)  is  found  from 


The  velocities  Uv(x),  uv(x,v,t),  equation  (23)  will  now  be 
written  in  a  more  useful  form  in  terms  of  the  free-stream  components 
V(x)  and  W(x),  and  the  radial  and  tangential  perturbation  velocities 
ur(r,e)  and  uQ(r,e). 

Uv(x,r,e)  =  W(x)  cos  (m-£)  +  V(x)  sin  (m-£)  (28) 

W(x),  V(x)  are  the  free-stream  component  contributions  in  the  v 
directions,  see  figure  4. 


uy(r,0)  =  ur(r,0)cos[0-(m  -  |)]  -  u0(r,0)sin[0-(m  -  J)](29) 


The  inner  potential  $  is  the  potential  that  satisfies  the  flow 

2 

tangency  condition  on  the  actual  body  surface.  This  assumes  that  the 

effects  of  the  other  two  terms  in  equation  (3)  cancel  one  another  on 

the  body  surface.  This  condition  holds  true  on  the  surface  of  the 

equivalent  body  and  is  satisfied  approximately  on  the  actual  body 

surface.  For  the  cases  studied  in  Ref  4,  the  inclusion  of  the  other 

two  terms  had  a  small  effect  (Ref  4:11).  The  effects  of  only  the  inner 

potential  <j>  will  be  considered  in  the  boundary  condition. 

2 


The  perturbation  velocities  in  equation  (29)  are  related  to  the 


crossflow  or  inner  potential  $  by: 

2 

ur  a<J>  (r,e) 

V  (r,e^  =  ~lr 

oo 

uQ(r, 8)  _  l  (r.e) 
V  r(e)  30 

Differentiation  of  equation  (11)  gives: 


MH  na„cos„e 
V  n  n 

n=l  rn+1(e) 


nansinne 

rn+1(0) 


( 


( 


( 


Substitute  equations  (31)  and  (32)  into  (29).  Using  this  result  and 
equation  (28),  the  flow  tangency  condition  given  by  equation  (23) 
becomes : 


The  perturbation  term  ux(x,v,t)  is  omitted  since  it  is  small 
compared  to  Ux(x).  The  rate  of  change  with  x  of  the  body  cross-section¬ 
al  area  S'(x)  can  be  written  in  terms  of  the  equivalent  body  radius  and 
thus  the  source  term  in  equation  (33)  can  be  written 


U  (x) 
2ttV 

00 


Ux(x) 

2i?V 

co 


£  (*Req>  - 


(34) 


where  Req  is  the  radius  of  the  equivalent  axisymmetric  body. 

Equation  (33)  is  rewritten  by  transposing  and  factoring  terms  and 
is  applied  at  MC  control  points.  The  result  is  a  set  of  MC  flow 
tangency  equations  in  MH  unknown  polar  harmonic  coefficients  an 


MH 

I 

n=l 


na. 


rn+1(6,.) 


-  cos  ne..  cos  [e..-(m..  -  |-)] 


(35) 


+  sin  ne.  sin  [e.  -  (m^  -  -|)]j 


Ux(x)  c  i  \ 

ST"  rfeT)  “«  -  in 

00  I 


Ux(x)  dv(0.) 
V  dx 


.  Mix) 


v  cos(mi- 
00 


i  *  1,2 . MC  ,  MC<MH 


The  problem  now  is  to  find  the  best  possible  values  for  a  set  of 
unknown  a^^,. .  .a^  for  a  set  of  MC  linear  equations.  Instead  of  the 
exact  solution  of  MH=MC,  a  set  of  MC  equations  will  be  satisfied  in 
the  least-square  sense.  Minimize  the  quantity 


Fi 


b 
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(36) 


N 


MC 

E=  I  6  2 

i=l  1 

where  6.. 2  represents  equation  (35).  (Equation  (35)  [A]=[B],  [A]  -  [B] 
=  <5..)  Using  the  procedure  (Ref  9:177-178)  to  minimize  E  by  setting 
the  partial  derivatives  equal  to  zero. 


and  letting 


!I  _SjL  8E 

8a  ’  8a  *  **’*  8aMU 

2  rlH 

*i  =  ei  -  mi +  i 


gives  the  result 


MH 

A 


MC  cos.e.cosit.  +  sin  .e.sinii). 
v  ji  i  Ji  i  _ n 

i=l  '  rj+1(9,)  rn+1 


(-cos  ne.  cos^  +  sin  ne^  sin  ^..) 


(37) 


(38) 


(39) 


Ux(x)  dv(e.) 
~V  dx~ 


W(x) 

V 


cos(m..  -  j ) 


Ux(x) 

■  2^r~ 
00 


n 


rn+l(8l) 


(-cos  nei  cos  ei  +  sin  ne^  sin  ) 


n=l,2, . . . ,MH 
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This  represents  a  set  of  MH  equations  in  inknowns  aj,  a2»...»aMH 
which  are  to  be  solved  simultaneously.  The  numbers  MC  and  MH  must  be 
specified,  and  their  determination  requires  judgement.  The  choosing  of 
MC  and  MH  will  be  discussed  in  detail  in  a  later  section,  in  which  the 
input  to  the  trajectory  program  is  determined.  The  details  are  in  (Ref 
4:15-18). 

Fuselage  Mounted  Air  Inlets  (Ref  4:18-20)  Air  inlets  include  the 
entire  region  from  the  mouth  to  the  engine  exhaust.  The  portion  of  the 

inlet  that  is  a  solid  boundary  is  treated  as  part  of  <f>  ,  the  inner 

2 

fuselage  potential. 

In  the  section  where  the  contour  is  not  a  solid  boundary,  figure  5, 
the  body  slope  boundary  condition  used  in  the  polar  harmonic  calculation 
is  determined  in  a  different  way.  The  inlet  is  first  considered  to  have 
a  solid  boundary.  Control  points  are  laid  out  on  the  range  0°<e<180°  of 

the  crossflow  contour.  With  each  control  poisiis  associated  a  stream- 

•  *  ^ 

wise  body  slope  and  upwash  if  angle  of  attack  is  considered.  For  the 
control  points  lying  on  the  contour  which  is  not  a  solid  boundary  the 
streamwise  slopes  0^  depend  on  the  inlet  to  free-stream  velocity  ratio, 
designated  Vg/V^.  For  Vg/V^l,  streamlines  are  parallel  to  the  axis  of 
the  inlet,  and  for  lesser  values  they  point  downward.  To  account  for 
the  change  in  streamline  direction  caused  by  blockage,  the  streamwise 
slopes  calculated  for  the  solid  boundary  to  give 


V  (1 


VD)  dv(e) 
V  ]  dx 

oo 


(40) 


solid  boundary 


The  term  in  the  flow  tangency  condition  in  equation  (33)  is 

replaced  by  e  .  The  rate  of  change  of  cross-sectional  area  S'(x)  must 

3 

also  be  determined.  For  a  Vg/V^O.5  only  one  half  of  the  inlet  cross- 
sectional  area  is  added  to  the  fuselage.  This  adjusts  the  area  dis¬ 
tribution  to  reflect  the  inlet  velocity  ratio. 

The  outer  potential  eg  is  also  adjusted  to  reflect  the  inlet 
velocity  ratio.  Its  cross-sectional  area  distribution  is  modified  by 
excluding  a  portion  of  the  inlet  area  based  on  Vg/V^. 

Wing-Pylon  Flow  Model  (Ref  4:20-24)  The  fuselage  is  first  repre¬ 
sented  by  the  potential  flow  methods  previously  described.  The  influ¬ 
ence  of  the  fuselage  on  the  exposed  wing  panels  is  then  determined. 

This  is  done  by  a  vortex  lattice  with  unknown  vortex  strengths  which  is 
laid  out  on  the  exposed  wing  panels  and  pylon.  An  image  vortex  lattice 
of  the  wing  is  constructed  inside  the  fuselage.  The  wing-pylon  loading 
is  computed  in  terms  of  the  vortex  strengths  with  the  inclusion  of  the 
fuselage  influence  on  the  exposed  wing  panels  and  pylon.  The  inclu¬ 
sion  of  the  image  vortices  inside  the  fuselage  approximately  satisfies 
the  condition  of  no  flow  through  the  surface  and  accounts  for  nearly 
all  of  the  wing-fuselage  interference.  The  crossflow  plane  theory  for 
the  imaging  scheme  is  described  in  (Ref  4:20-23)  and  (Ref  9). 

The  vortex  lattice  is  set  up  only  on  the  portion  of  the  wing  out¬ 
side  the  equivalent  body  of  revolution  for  the  noncircular  fuselage. 

If  t.he  body  of  revolution  radius  varies  over  the  region  of  influence 
of  the  wing,  the  maximum  radius  is  used.  Each  horseshoe  vortex  that 


is  outside  the  body  of  revolution  is  then  imaged  inside.  Velocities 
normal  to  the  equivalent  body  are  exactly  zero  in  the  wing  chordal  plane 
and  are  very  small  at  all  other  locations  on  the  body  of  revolution  sur¬ 
face  near  the  wing-body  junction.  For  noncircular  bodies  the  actual 
fuselage  surface  is  generally  close  to  the  equivalent  body  surface  so 
that  to  a  good  approximation  the  same  conditions  hold  for  the  actual 
fuselage  surface. 

Vortex  Lattice  (Ref  4:20-25)  Figure  6  shows  the  vortex  lattice 
arrangement.  The  wing  panels  and  pylons  are  divided  into  trapezoidal 
area  elements.  A  horseshoe  vortex  is  placed  in  each  area  element  such 
that  the  spanwise  bound  leg  lies  along  the  element  quarter  chord  and 
its  trailing  legs  along  the  sides  of  the  element.  The  trailing  legs 
lie  in  the  plane  of  the  area  element.  The  area  elements  in  each  chord- 
wise  row  have  equal  chords  and  spans.  In  the  spanwise  direction,  the 
area  element  width  need  not  be  equal  to  allow  for  closer  spacing  where 
large  spanwise  loading  gradients  exist.  For  a  wing  with  breaks  in 
sweep  and/or  dihedral,  the  area  elements  are  arranged  spanwise  so  that 
the  breaks  lie  on  the  line  formed  by  one  of  the  sides  of  a  chordwise  row 
of  elements.  The  wing-pylon  junction  is  also  made  to  lie  along  a  common 
boundary  between  two  adjacent  rows  of  elements.  On  the  wing,  the  flow 
tangency  boundary  condition  is  applied  at  a  set  of  control  points  given 
by  the  midpoint  of  the  3/4  chord  line  of  each  area  element  located  in 
the  wing  chordal  plane.  This  is  the  planar  approximation.  On  the 
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camberless  pylon,  the  flow  tangency  condition  is  applied  at  control 
points  situated  in  its  chordal  plane  determined  by  halving  the  thick¬ 
ness  envelope. 

Figure  7  illustrates  the  boundary  condition.  The  velocities 
normal  to  the  wing  consist  of  a  component  of  the  free-stream,  pertur¬ 
bation  velocities  u,  v,  and  w  induced  by  the  wing-pylon  horseshoe 
vortex  system,  and  perturbation  velocities  u..,  v^ ,  and  w^  induced  by 
the  distribution  of  sources  accounting  for  thickness  superimposed  on 
the  vortex  lattice.  This  thickness  distribution  will  be  well  described 
in  the  next  section.  With  M  control  points  on  the  left  wing  panel  and 
MP  on  the  left  pylon,  the  boundary  condition  on  the  left  wing  panel 
is  given  by 
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(41) 


v  =  1,2, ... ,M 

With  the  camberless  pylon  at  zero  incidence,  the  pylon  boundary 
condition  is  written  for  MP  control  points 
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V  =  M+l,  M+2 . M+MP 

The  wing  angle  of  attack  is  a  and  a is  the  local  angle  of  attack 
due  to  wing  camber  and  twist.  Both  angles  are  assumed  small. 

The  right-hand  side  of  equation  (41)  represents  the  free-stream 
component  and  the  externally  induced  perturbation  velocities  normal  to 
the  wing  chordal  plane.  The  right-hand  side  of  equation  (42)  consists 
only  of  an  external  perturbation  velocity  normal  to  the  pylon  chordal 
plane  since  there  is  no  pylon  incidence.  The  first  summation  on  the 
left-hand  side  of  the  equations  represents  the  perturbation  velocities 
induced  by  vorticity  on  the  left  and  right  wing  panels.  The  second 
summation  represents  the  velocities  induced  by  vorticity  on  the  left 
and  right  pylons. 

The  functions  Fu>  Fy,  Fw,  are  influence  coefficients  relating  the 
perturbation  velocity  components,  induced  at  some  point  by  a  horseshoe 
vortex,  to  its  circulation,  and  the  coordinates  of  the  point  relative 
to  the  vortex.  Functions  F.  and  F.  relate  the  perturbation  velocity 
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components,  induced  at  some  point  by  a  wing  image  horseshoe  vortex,  to 
its  circulation  and  the  coordinates  of  the  point  relative  to  the 
bound-leg  midpoint  of  the  image  horseshoe  vortex. 


Using  the  Biot-Savart  Law  (Ref  11:156-160)  the  perturbation 
velocities  induced  at  a  point  by  a  horseshoe  and  image  horseshoe  vortex 
are 


u(x,y,z)  =  5ST  CFu(x.y,z,S,*,+)  +  Fju(x,y,z,sjv,*jv,<,1.v)] 


x(x.y.x)  ■  4^r  [Fv(x,y,z,s,*,*)  +  Fjv(x,y,z,s.  ,*,,*.  )] 

ao 

w(x,y,z)  =  jjj-  [Fw(x,y,z,s, ![/,<(.)  +  Fiw(x.y •z,siv,4»iv»<J>iv)]  (43) 


The  influence  functions  depend  only  upon  the  coordinates  (x,y,z) 
of  the  point  at  which  the  velocity  is  to  be  computed  relative  to  the 
horseshoe  vortex/image  horseshoe  vortex,  the  vortex  span,  s,  the  image 
vortex  semispan  s^,  bound-leg  sweep  angle  ^iu,  and  dihedral  angle 
<j> , i u •  The  analytical  expressions  for  the  influence  coefficients  are 
given  (Ref  2:12-14),  and  are  subject  to  equation  (13),  (Ref  2:14)  to 
account  for  the  images  of  the  horseshoe  vortices  laid  out  on  the  right 
wing  panel . 

Equations  (41)  and  (42)  represent  a  set  of  M  +  MP  simultaneous 
equations  in  which  the  unknowns  are  the  M  +  MP  values  of  circulation 
strength  r.  These  can  be  solved  for  a  given  angle  of  attack,  twist  or 
camber  distribution,  and  a  specified  set  of  externally  induced  pertur¬ 
bation  velocities  u./V^,  v^/V^,  and  w^V^  caused  by  wing  and  pylon 


thickness  effects  and  from  other  aircraft  components  such  as  fuselage, 
rack  and  stores. 

Thickness  Model,  Ming  and  Pylon  (Ref  4:25)  Wing  and  pylon  thickness 
are  accounted  for  by  laying  out  thickness  strips  on  the  exposed  panels. 
These  strips  are  made  of  a  three-dimensional  source  distribution. 

The  three-dimensional  source  method  is  similar  to  that  used  by 
Fernandes  (Ref  12).  Figure  8  and  9  show  the  coordinate  and  angle  sys¬ 
tem  for  the  wing  and  pylon. 

The  incompressible  velocity  potential  due  to  a  surface  of  sources, 
with  chord  |Xb-Xa|  and  span  |Yb-Ya|,  on  the  wing  may  be  obtained  (Ref  6; 
Ref  13). 


M  -  _L 

V  2 it 

oo 


tan  <j>t  dXj  dYj 

C(X1-X)2+  (YrY)2  +  (ZrZ)2]* 


(44) 


for  a  surface  located  at  Z^constant  and  local  incidence  angle  <j>t 
measured  from  the  positive  x  direction.  The  integration  is  performed 
over  the  span  and  chord  of  the  surface  and  coordinates  X,  Y,  Z  locate 
the  field  point  at  which  the  velocity  potential  is  determined.  The 
integration  is  simplified  and  performed.  This  gives  A^/V^  for  one 
strip.  This  is  differentiated  to  give  the  perturbation  velocities  due 
to  one  thickness  strip.  The  same  method  is  applied  to  the  pylon  (Ref  2: 
15-20) 
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The  perturbation  velocities  induced  by  the  thickness  envelopes  of 
the  wing  and  pylons  are  computed  at  the  control  points  associated  with 
the  vortex  lattice  area  elements  on  the  wing-pylon  combination.  The 
velocities  are  obtained  by  adding  the  result  of  summing  equations  over 
all  wing  source  strips  to  the  result  of  summing  equations  over  all  pylon 
source  strips.  The  final  sums  are  then  added  to  the  set  of  externally 
induced  perturbation  velocities  that  appear  in  the  flow  tangency  bound¬ 
ary  condition.  All  interference  effects  induced  by  the  wing  on  the 
pylon  and  vice-versa  are  accounted  for  and  will  be  reflected  in  the 
resulting  circulation  strength  distribution  on  the  wing-pylon  combina¬ 
tion  (Ref  2:15-20). 

Store  and  Rack  Flow  Models  (Ref  4:26)  The  store  model  uses  the 
same  technique  as  the  axisymmetric  fuselage  model.  The  model  accounts 
only  for  the  volume  distribution  of  the  axisymmetric  body,  by  a  dis¬ 
tribution  of  three-dimensional  point  sources  along  the  longitudinal 
axis  of  the  store. 

The  rack  model  is  approximated  by  an  axisymmetric  body  with  three 
short  pylons  attached.  The  body  is  modeled  by  a  source  distribution 
while  the  short  pylon  is  neglected  (Ref  2:20). 

Force  and  Moment  Calculation  (Ref  4:26-27) 

Calculation  of  the  forces  and  moments  on  a  store  requires  speci¬ 
fication  of  the  nonuniform  velocity  field  in  which  the  store  is 
immersed.  This  field,  as  seen  with  respect  to  the  store,  must  be 


determined  at  each  point  in  time.  The  field  includes  the  free-stream 
velocity,  the  perturbation  velocities  induced  by  the  parent  aircraft 
and  the  angular  velocities  due  to  the  store's  pitch,  yaw,  and  roll 
motions.  The  velocity  distribution  along  the  body  axis  is  required 
to  calculate  the  body  forces,  and  the  velocity  distribution  over  the 
tail  fin  surfaces  is  required  for  the  empennage  forces  and  moments. 


Velocity  Field  (Ref  2:21-26) 


Nomenclature: 


P«q»r  -  rotational  velocities  about  x,y,z  axes 

r  -  radial  distance  in  y-z  plane 

us’vs’ws  "  sum  perturbation  velocities 

U°oS  x,V<*>c  x/»W«»c  7“  free-stream  velocity  components  as  seen 
’  s,y  s,z  by  store>  f.gure  1Q 

Us,Vs,Ws  -  total  velocities  as  seen  by  store 

*SjItl  -  x  location  of  store  moment  center 

5>n,c  -  inertial  coordinate  system  fixed  in 

fuselage  nose,  positive  forward  along 
longitudinal  axis,  positive  laterally 
to  the  right,  and  positive  vertically 
downward 


The  x  ,  ys,  zs  coordinate  system  to  be  used  is  shown  in  figure  10. 
The  coordinate  system  is  fixed  in  the  store,  with  the  origin  at  the 
store  nose.  The  velocities  Us,Vg,Ws  are  total  velocities  as  seen  by  a 
point  on  the  store  and  each  is  composed  of  three  components:  free- 
stream,  perturbation  due  to  parent  aircraft,  and  damping  due  to  store 


motion. 
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A  reference  system  must  now  be  established  to  represent  the 
equations  of  motion  used  in  the  trajectory  calculation.  Figure  11 
shows  the  £,n,s  system  which  is  an  inertial  system  fixed  in  the  air¬ 
craft  which  is  in  rectilinear  flight  at  uniform  velocity.  The  positive 
£  axis  is  forward  along  the  fuselage  longitudinal  axis,  the  n  axis  is 
positive  to  the  right,  and  the  ?  axis  is  positive  downward.  The 
origin  is  fixed  in  the  nose  of  the  aircraft. 

In  figure  11,  the  origins  of  the  two  systems  have  been  drawn  to 
coincide  in  order  to  show  the  angles  used  to  determine  the  orientation 
of  the  store  with  respect  to  the  inertial  (£,n,0  axis.  The  Euler 
angles  relate  the  two  coordinate  systems  (Ref  14:100-103). 
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with  A  given  by  equation  (28) (Ref  2:22). 

For  a  parent  aircraft  moving  at  velocity  V^,  and  flying  at  angle 
of  attack  c^,  the  velocity  of  the  inertial  coordinate  system  relative 
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to  a  point  in  space  is 

y»=  v«COSafe?  +  ^sinafe?  (47 

The  velocity  of  the  store  moment  center  relative  to  the  moving  inertial 
system  is 

7  =  le.  +  ne  +  ££  (48 

m  £  n  c 

The  velocity  of  the  store  relative  to  a  point  in  space  is 

V’- ♦  ’«,  <49 

Since  the  velocity  components  needed  in  equation  (45)  are  with 
respect  to  a  point  fixed  in  the  store  and  moving  with  the  store, 
equation  (49)  becomes 

VccS  =  -  (V^cosc^  +1)5^  -  nin  (50 

-  (Vaosinaf  +i)e{. 

The  components  in  the  x,y,z  coordinate  system  of  figure  10  are  then 


f“  - 

U”s,x 

'  -  (VooCOSaf  +|) 

v»s,y 

-  [A]T 

• 

-  n 

.Woos,z . 

(V.sinaf  +c). 

where 


U“s,xs  =  "U"s,x  (52) 

V  =  V 
“s,ys  “s,y 

W  =  -U 

“s,z$  °°s,z 


It  is  essential  to  remember  that  the  second  term  in  each  equation, 
equation  (45),  is  comprised  of  the  perturbation  velocities  induced  by 
the  fuselage,  wing,  pylon,  rack,  and  other  stores.  These  velocities, 
determined  in  the  previous  sections,  apply  to  the  equivalent  incom¬ 
pressible  configuration.  At  each  point  in  the  trajectory  the  points 
at  which  the  velocities  are  required  must  be  located  in  the  incompres¬ 
sible  space.  The  velocities  calculated  using  specific  equations  for 
each  component  are  defined  in  the  coordinate  system  of  each  component. 
Prior  to  transforming  these  back  to  the  compressible  space  they  must  be 
summed  in  the  fuselage  coordinate  system  since  this  is  the  coordinate 
system  in  which  the  compressibility  correction  was  applied  (Ref  4:26- 
27).  This  is  the  U.n.t)  system  of  Figure  11. 

Let  ul,v',w'  be  the  sums  in  the  £»n,c  directions.  From  the 
%  n  t 

section  on  compressibility  correction 
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(53) 
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The  velocity  components  in  the  x,y,z  direction  of  figure  10  are 


u 

- 

1  u. 
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Cj 

(54) 


and  finally,  the  components  in  the  coordinate  system  of  figure  10  are 


u$  =  -u  ,  vs  =  v  ,  w$ 


=  -w 


(55) 


The  velocities  to  be  used  in  the  force  and  moment  calculations 
are  nondimensionalized  by  the  store  free-stream  velocity.  Thus,  equa¬ 
tion  (45)  becomes 


(56) 
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where  from  equation  (49)  or  (50) 


Vo»s  =  [(Voocosaf  +  l)2  +  n2  +  (Voosinaf  +£)2]*  (57) 

At  any  point  in  time  during  the  trajectory,  the  forces  and  moments 


are  determined  by  removing  the  store  from  the  flow  field,  determining 


>  J  s 


the  velocities  at  a  series  of  points  which  the  store  longitudinal  axis 
and  tail  fins  occupied,  and  then  immersing  the  store  in  this  flow  field. 

Buoyancy  Forces  and  Moments  (Ref  2:27-29)  The  buoyancy  forces  and 
moments  are  calculated  with  the  assumption  that  the  flow  in  planes 
perpendicular  to  the  axis  of  the  store  obeys  Laplace's  equation  and 
that  the  upwash  and  sidewash  velocities  vary  along  the  store  length. 

The  potential  is  constructed  on  this  basis,  then  integration  of  the 
body  pressures  obtained  from  the  unsteady  Bernoulli  equation  yield 
both  the  buoyancy  forces  and  the  slender-body  forces.  The  final 
expression  for  the  buoyancy  normal-force  coefficient  is  given  by 
equation  (63)  of  reference  1.  The  derivation  is  carried  out  in  Ap¬ 
pendix  I  of  that  reference  with  the  final  expression  which  contains 
both  buoyancy  and  slender-body  terms  given  by  equation  (1-14). 

The  expression  used  to  calculate  the  coefficients  are 
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The  expressions  for  pitching  moment  and  yawing-moment  coefficients 
with  the  moment  taken  about  x$  m  (see  figure  10),  are 
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Slender-Body  Forces  and  Moments  (Ref  2:28)  The  slender-body  forces 
and  moments  are  derived  in  section  6.21  of  reference  1.  The  expressions 


(cN)  =  r1  rs*°  df  (aV)dx< 
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The  integrals  above  are  assumed  valid  until  the  separation  loca¬ 
tion  is  reached.  At  that  point,  the  viscous  forces  become  important. 

There  are  no  definite  rules  for  picking  x$  Q,  the  upper  limit  of 
integration  where  the  integrals  are  valid.  However,  the  problem  is 
discussed  In  section  6.2,  and  8.2  of  reference  1.  From  that  reference, 
a  formula  is  given  which  relates  x  /£  ,  which  is  the  location  on  the 

«  y  U  J 

body  where  potential  flow  is  no  longer  applicable,  to  x  ./£ 

J  )  1  b 

X  X 

=  0.378  +0.527  (66) 

s  s 


where  xg  ^/£s  is  the  location  on  the  body  of  the  maximum  negative  rate 

of  change  of  cross-sectional  area.  This  can  be  used  to  estimate  the 

value  of  x,  . 

s,o 

Viscous  Crossflow  Forces  and  Moments  (Ref  2:29-30)  From  the  as¬ 
sumed  separation  location  to  the  base  of  the  store  a  viscous  crossflow 
calculation  is  used  in  place  of  the  slender  body  calculation.  The 
expressions  are  derived  in  section  6.2.2  of  reference  1. 
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where  \IQ*  =  [(Vg*)2  +  (W$*)2]i 
r  drag  per  unit  length 

Ldc  ‘  qoos2a 
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Empennage  Forces  and  Moments  A  method  of  calculating  the  forces 


and  moments  associated  with  planar  or  cruciform  empennages  is  presented 
in  section  5.3,  reference  2.  The  details  are  not  essential  to  the 
construction  of  the  trajectory  program  input  and  will  not  be  presented 
here. 

Equations  of  Motion  (Ref  4:27;  Ref  2:33-38) 

The  complete  derivation  of  the  equations  of  motion  of  a  rigid  body 
with  mass  and  inertia  asymmetries  is  presented  in  Appendix  II  (Ref  2:131 
145).  Generally,  a  store  body  will  have  axes  of  geometrical  symmetry 
about  which  the  store  forces  and  moments  are  determined.  By  mass  asym¬ 
metry  it  is  meant  that  the  store  center  of  mass  does  not  have  to  lie  at 
the  origin  of  this  coordinate  system.  By  inertia  asymmetry  it  is 
meant  that  the  principal  axes  of  inertia  of  the  store  do  not  coincide 
with  the  geometric  axes  of  symmetry  so  that  the  tensor  of  inertia 
possesses  off  diagonal  terms,  the  products  of  inertia. 

The  derivation  of  the  equations  of  motion  used  two  previously  dis¬ 
cussed  coordinate  systems,  see  figure  11.  The  x,y,z  coordinate  system 
is  fixed  in  the  store  and  rotates  with  the  store.  The  x-axis  is  posi¬ 
tive  forward,  y-axis  positive  to  the  right,  and  the  z-axis  positive 
downward.  These  axes  coincide  with  the  geometric  axis  of  symmetry  with 
the  origin  at  the  point  about  which  the  moments  and  products  of  inertia 
and  the  aerodynamic  moments  are  calculated.  The  x,y,z  coordinate  system 


is  used  because  the  time  derivatives  of  the  moments  and  products  of 
inertia  do  not  appear  in  the  rotational  equations  of  motion.  The 
equations  of  motion  are  given  in  Appendix  II  (Ref  2:136,  144)  by  equa¬ 
tions  (11-15)  and  (11-40). 

The  x,y,z  coordinate  system  does  not  allow  the  store  position 
and  orientation  to  be  determined  relative  to  the  host  aircraft.  The 
C,n,e  system  is  used  for  this  purpose,  with  the  assumption  that  the 
parent  aircraft  is  flying  at  constant  velocity,  constant  angle  of 
attack,  and  constant  flight  path  angle  relative  to  the  horizontal.  The 
C»n,s  system  can  be  regarded  as  inertial  because  it  is  non-rotating, 
and  the  store  motion  relative  to  the  moving  system  calculated. 

The  5,n »C  system  is  fixed  at  the  nose  of  the  parent  aircraft  with 
the  c-axis  positive  forward,  n-axis  positive  to  the  right,  and  c  posi¬ 
tive  vertically  downward. 

The  orientation  of  the  store  with  respect  to  the  £,n,c  system  is 
accomplished  through  Euler  angles  v,0,4  (Ref  14:100-107).  The  time 
histories  of  these  angles  are  determined  by  expressing  them  as  store 
rotational  velocities  p,q,r,  see  figure  11.  The  differential  equations 
expressing  this  are 

$  =  (qsin$  -  rcos4>)/cos0 

0  =  (qcos$  -  rsin$)  (73) 


$  =  p  +  qsin$  tano  +  rcos$  tan© 


These  equations  along  with  equations  (11-16)  through  (11-18),  and  (II- 
41)  through  (11-43)  are  the  final  differential  equations  which  are  to 
be  integrated  to  determine  the  store  position  and  orientation  relative 
to  the  inertial  system  fixed  in  the  fuselage  (Ref  2:136-144). 


In  the  equations,  S.n,?  are  the  accelerations  of  the  store  moment 
center  in  the  inertial  coordinate  system.  The  velocities  in  this  sys¬ 
tem  are  |,n,c  and  the  quantities  xQ,  yQ,  zQ  are  these  velocities  re¬ 
solved  into  the  store  body  coordinate  system  using  [A]^,  previously 
defined.  The  location  of  the  store  center  of  mass  relative  to  the 
store  moment  center,  in  store  body  coordinates,  is  x,  y,  z.  The  rota¬ 
tional  velocities  in  the  store  body  coordinate  system  are  p,  q,  r. 

The  quantities  Ivv,  I,  I,,  are  the  moments  of  inertia  and  I  ,  I  , 
xx  yy  zz  yz  xz 

I  are  the  products  of  inertia, 
xy 

The  remaining  parameters  in  the  equations  of  motion  are  the  store 
mass,  m,  the  forces  Fx,  Fy,  Fz;  the  moments,  Mx,  My,  Mz.  The  forces 
and  moments  act  on  the  store  and  are  positive  in  the  positive  x,y,z,p, 
q,r  directions. 

Fx  *  m9x  -  9-s  SR  CA 


Fy  =  m9y  +  9a>s  Cy 


Fz  =  mgz  "  q»s  SR  CN 


where  g  ,  g  ,  g  are  the  components  of  the  gravitational  acceleration 
a  y  z 

in  the  store  body  coordinate  system  (Ref  2:36).  The  second  terms  in 
equation  (74)  are  the  aerodynamic  forces  acting  on  the  store.  The 
axial-force  coefficient,  C^,  is  specified  by  the  user  and  is  positive 
in  the  negative  x  direction.  The  sideforce  and  normal  force  coefficients, 
Cy  and  are  calculated  by  the  computer  program  using  the  methods  de¬ 
veloped  in  the  section  on  force  and  moment  calculation.  The  parameters 
q<x.s  and  SR  are  the  dynamic  pressure  and  reference  area  used  in  nondimen- 
sionalizing  the  forces. 


q»s  =  ip»V 

SR  =  ^max 


(75) 

(76) 


The  density  is  assumed  constant  at  the  value  for  the  aircraft 
flight  altitude  at  time  t=0,  and  Va>s  is  given  by  equation  (56).  The 
maximum  store  body  radius  is  used  for  a 

max 

The  moments  about  the  store  moment  center  are 


MX =  <K  Vrc*  +  m<v  -  y> 

My  *  '’“s  Win  +  '  9Z;) 

Hz  *  q-s  SR*RCn  +  "(Sy5  *  9X*> 
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The  coefficients  C  ,  C  ,  C  are  calculated  by  the  computer.  The 
X*  HI  n 

reference  length  is  the  maximum  store  diameter 


(78) 


The  second  terms  in  the  moment  equations  are  the  moments  produced  by 
the  gravitational  force  when  the  location  of  the  store  center  of  mass 
does  not  coincide  with  the  center  of  moments. 

The  final  specifications  required  before  the  equations  of  motion 
can  be  integrated  are  the  initial  conditions.  These  are  the  position 
of  the  store  center  of  moments  (c,n,s),  the  translational  velocity  of 
this  point  (£,n,£)»  the  angular  velocity  about  the  store  axes  (p,q,r) 
and  initial  orientation  angles  s',©,#.  With  these  specifications  the 
equations  can  be  integrated  with  the  aerodynamic  forces  and  moments 
being  recalculated  at  each  point  in  the  trajectory. 

The  mathematical  model  for  determining  the  flow  field  around  an 
aircraft  fuselage-wing  combination  and  associated  stores  has  been 
developed.  This  model  has  been  used  to  establish  the  forces  and  moments 
on  a  separated  store,  and  using  these  forces,  the  equations  of  motion 
can  now  be  integrated  to  predict  the  store  trajectory. 

The  equations  and  relationships  have  been  programmed  into  a  source 
and  trajectory  program  (Ref  5).  The  input  for  each  of  these  programs 
will  now  be  described  and  developed. 
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Figure  7.  Flow  Tangency  Boundary  Conditions 
(Ref  2:  61) 
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Figure  8.  Coordinate  System  for  Wing  Thickness 
(Ref  2:  62) 


PYLON  THICKNESS 
ENVELOPE 


Figure  9.  Coordinate  System  for  Pylon  Thickness  (Ref  2:  63) 


Figure  10.  Coordinate  System  Fixed  in  Ejected  Store 

(Ref  4:  64) 
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III.  NE ILSEN  PROGRAM  APPROACH 


The  Neil  sen  approach  uses  two  separate  computer  programs,  a  source 
program  and  the  trajectory  program  (Ref  5:4-8).  The  two  programs  are 
used  with  no  direct  exchange  of  information  during  a  run.  The  source 
program  does  provide  an  input  to  the  trajectory  program,  but  the  user 
must  provide  the  information  exchange. 

The  source  program  is  used  to  represent  an  axi symmetric  body  as  a 
distribution  of  sources  along  the  axis  of  the  body.  This  representation 
is  used  for  the  fuselage,  store,  and  pylon  volume  effects  on  the  overall 
flow  field  around  the  aircraft.  The  program  calculates  and  prints  the 
source  strengths  and  locations.  These  quantities  are  used  as  input  data 
to  the  trajectory  program. 

The  trajectory  program  calculates  the  trajectory  of  the  dropped 
store  by  accounting  for  the  effects  of  the  other  aircraft  components  on 
the  store,  and  integrating  the  equations  of  motion.  The  program  con¬ 
tains  options  to  allow  for  a  circular  or  noncircular  fuselage,  no 
fuselage,  pylons,  racks,  and  up  to  ten  stores.  This  report  will  deal 
only  with  a  noncircular  cross-section  fuselage,  with  one  store  on  the 
fuselage  and  one  store  separated.  There  is  no  rack  or  pylon.  This 
configuration  was  chosen  to  compare  with  available  wind  tunnel  data. 

The  input  for  each  program  will  be  discussed  in  detail  in  Appendix 
B  and  Appendix  C.  Those  items  which  do  not  apply  will  be  referenced  to 
the  proper  users  manual.  Details  of  the  preparation  are  included  to 
allow  future  users  to  concentrate  on  those  items  which  cause  the  most 
difficulty. 


IV.  Source  Program  Input  and  Analysis 


The  program  input  (Ref  5:9-16)  is  initiated  by  the  user  represent¬ 
ing  the  body  as  an  equivalent  body  of  revolution  (EBR).  The  surface  of 
the  EBR  is  then  approximated  by  a  series  of  polynomials  which  represent 
the  EBR  x,r  distribution.  At  this  point  the  source  program  is  run  with 
a  user  specified  finite  source  distribution.  The  program  calculates 
the  shape  of  the  EBR  from  the  source  strengths  and  the  output  is 
presented  so  that  the  user  may  compare  calculated  values  of  the  surface 
to  the  polynomial  representation  of  the  surface.  There  is  no  set  pro¬ 
cedure  for  specifying  the  source  distribution;  trial  and  error,  and 
experience  are  the  only  tools  available.  Once  the  proper  representa¬ 
tion  is  determined,  the  source  strengths  and  locations  are  input  into 
the  trajectory  program. 

The  source  program  input  was  developed  for  the  fuselage  and  store 
at  M  =  0.6  and  M  =  0.9.  The  fuselage  input  will  be  described  in  detail. 
The  input  was  determined  in  distinct  steps: 

1.  Cross-section  area  distribution  for  fuselage 

2.  Inlet  velocity  ratio 

3.  Cross-section  area  distribution  of  EBR 

4.  Polynomial  to  represent  EBR 

5.  Develop  input  deck 

6.  Run  source  program 


Cross-Section  Area 


Cross-section  data  were  available  at  22  fuselage  stations.  Data 
on  more  stations  would  have  improved  the  program  accuracy  in  represent¬ 
ing  the  body.  The  area  at  each  station  is  converted  to  an  EBR  area  and 
radius.  This  provides  one  point  on  a  curve  which  must  then  be  repre¬ 
sented  by  a  given  polynomial.  The  accuracy  of  the  polynomial  is 
directly  related  to  the  number  of  data  points. 

Inlet  Velocity  Ratio 

As  previously  stated,  the  inlet  velocity  ratio  is  used  to  adjust 
the  area  distribution  to  reflect  the  change  in  streamwise  flow  caused 
by  blockage  due  to  the  presence  of  the  inlet.  For  an  inlet  velocity 
ratio  of  0.5  only  one-half  of  the  inlet  cross-sectional  area  is  added 
to  the  area  of  the  body.  For  an  inlet  velocity  ratio  of  1.0,  none  of 
the  cross-sectional  area  is  added.  This  correction  is  made  on  the 
right  and  left  nacelle. 

The  inlet  ratio  was  determined  from  experimental  data  (Ref  15)  to 
be  V^V.  =  0.862.  This  was  calculated  by  knowing  the  theoretical  cap¬ 
ture  area  of  the  duct,  the  exit  area  of  the  nozzle,  and  experimental 
flow  data.  By  assuming  a  constant  mass  flow  through  the  duct,  the 
velocity  at  the  duct  inlet  was  determined.  This  duct  inlet  velocity 
divided  by  the  free-stream  velocity  gives  the  required  inlet  velocity 
ratio.  Reference  5  uses  an  inlet  velocity  ratio  of  0.5  with  no  ex¬ 
planation  of  how  the  number  was  determined.  This  value  could  be  used 
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if  no  experimental  data  were  available.  Note:  the  inlet  velocity 
ratio  must  be  greater  than  or  equal  to  zero,  but  less  than  or  equal  to 
1.0. 

O 

A  capture  area,  including  both  nacelles,  of  1068  in  was  used. 

With  an  inlet  velocity  ratio  of  0.862,  the  area  at  each  fuselage  sta- 

2 

tion  was  reduced  by  921  in'  . 

Cross-Sectional  Area  Distribution  of  EBR 

The  cross-sectional  area,  corrected  for  inlet  velocity  ratio  was 
determined  for  each  fuselage  station.  This  area  was  then  converted  to 
an  EBR  radius.  Table  2,  Appendix  A. 


EBR 


total 


1/2 


(79) 


The  information  on  the  last  station  was  taken  from  the  engineering 
drawing. 

Polynomial  to  Represent  EBR 

The  EBR  is  represented  by  a  series  of  polynomials  describing  the 
curve  of  r/i  plotted  versus  x/i,  figure  12.  All  coordinates  are  non- 
dimensional  ized  with  respect  to  the  body  length. 

The  equation  used  to  describe  the  EBR  is  given  by 

r/l  =  C  +  C  /C  (x/l)2  +  C  (x/£)2+  c  +  C  [x/l)**  C  [x/i)2  (80) 

1  7  2  3  4  5  6 


Combinations  of  coefficients  thru  C?  are  used  to  describe  the 
desired  curve.  For  this  case,  let  r/x,  =  y  and  x/£  =  x.  The  following 
curves  can  be  represented  by 

*  3  C!  +  V  C2X2  +  C3X  +  C4  +  C5X  +  C6x2  (81) 

Straight  line 

y  =  through  0=0 

y  ■  c  +  c  x  c=c=c=c=c=o 

1  5  2  3  i*  6  7 

Circle 


y  =  c  +  c  /  c  x z  +  c  x  +  c  c  =  c  =  o 

1  7  2  3  4  6  5 

Parabol  a 

y  -  C  +  C  x  +  c  x8  C  =C  =C  =  C  =0 

1  5  6  2  3  4  7 

Ellipse 

y  =  C  /  C  x2  +  C  C  =  C  =  C  =  C  =0 

7  2  4  13  5  6 

C  <  0 
2 

Hyperbola 

y  =  C  /  C  x2  +  C  C  =C  =  C  =  C  =0 

7  2  **  1  3  5  6 

C  <  0 


54 


The  source  program  allows  for  seven  polynomials  (NSECT)  to  be  used 
to  describe  an  EBR.  This  includes  a  polynomial  to  close  the  body,  or 
model  the  wake.  The  wake  was  modeled  (Ref  5; 13, 16)  for  the  ogive  store 
with  the  same  polynomial  as  the  nose.  For  the  present  case  a  straight 
line  was  used  to  close  the  fuselage  body  (EBR). 

Observation  of  the  EBR  plot,  figure  12,  indicated  that  straight 
lines  or  parabolas  would  be  sufficient  to  model  the  EBR.  A  least  square 
curve  fit  was  used  to  determine  the  coefficients.  Table  3,  Appendix  A. 
Source  Program  Input 

The  development  of  the  source  program  input  is  given  in  Appendix  B. 
The  two  variables  in  the  input  deck  are  NRAT  and  PERCR.  NRAT  is  the 
number  of  segments  which  the  body  will  be  divided  for  the  specification 
of  the  source  distribution.  PERCR  is  the  source  spacing  for  each  NRAT 
segment.  PERCR  is  input  as  a  fraction  of  the  local  body  radius  of  the 
segment. 

Source  Program  Analysis 

The  output  of  the  source  program  is  given  in  figure  29,  Appendix  B. 
The  data  from  this  program  were  input  into  the  trajectory  program. 

The  value  of  NRAT  was  chosen,  see  figure  12,  which  left  the  values 


of  PERCR  as  the  only  variables  in  the  source  program.  It  is  recommended 
(Ref  5:13)  that  an  initial  run  be  made  with  PERCR  =  1.0.  The  value  of 
PERCR  should  then  be  decreased  until  the  proper  shape  is  calculated  or 


a  maximum  number  of  100  sources  is  reached.  Computation  in  the  tra¬ 
jectory  program  is  a  function  of  the  number  sources  which  must  be  used 
to  model  the  fuselage,  stores,  and  racks.  The  maximum  of  100  sources 
is  applied  to  each  component  modeled. 

The  accuracy  of  the  program  generated  shape  was  judged  by  analysis 
of  the  output  section  (figure  29):  SHAPE  CALCULATED  FROM  SOURCE  DIS¬ 
TRIBUTION  AND  POLYNOMIALS.  For  each  x/2,  specified,  a  value  of  R/L 
(S.D.)  and  R/L  (POLY)  is  calculated.  A  guideline  was  established  for 
judging  the  acceptable  error  between  R/L  (S.D.)  and  R/L  (POLY).  Except 
for  the  nose  and  tail,  the  error  must  be  less  than  two  percent. 

Runs  were  made  with  PERCR  equal  to  1.0,  0.8,  0.6,  0.4.  The  value 
of  0.4  gave  more  than  100  sources.  Results  were  outside  the  guidelines 
for  all  of  these  values.  By  reference  to  figure  12,  it  was  noted  that 
the  body  radius  in  segments  one  and  two  was  small  compared  to  the 
others.  Also  the  accuracy  did  not  vary  with  values  of  PERCR.  PERCR 
was  held  at  1.0  in  these  segments  while  the  others  were  varied.  PERCR 
was  reduced  in  a  random  manner  in  segments  three,  four,  and  five  until 
values  of  0.9,  0.3,  and  0.4  respectively  were  found  to  give  the  best 
results. 

The  computer  execution  time  for  the  fuselage  source  program  with 
73  sources  was  0.96  seconds.  For  the  ogive  store  with  59  sources,  the 
execution  time  was  0.779  seconds.  These  values  were  taken  from  runs 
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V.  TRAJECTORY  PROGRAM  -  INPUT  AND  ANALYSIS 
Selection  of  the  Number  of  Polar  Harmonic  Coefficients 

The  trajectory  program  has  numerous  parameters  that  must  be 
either  calculated  or  set  before  the  program  can  be  run.  Most  para¬ 
meters  are  set  by  the  flight  condition  or  geometry  of  the  complete 
vehicle.  The  only  parameter  which  must  be  picked  on  the  basis  of  the 
trajectory  program  output  is  the  number  of  polar  harmonic  coefficients 
input,  item  13,  (see  Appendix  C).  Other  parameters  affect  the  output, 
but  they  are  set  by  the  given  problem.  The  polar  harmonics  must  be 
chosen  by  the  user,  and  very  few  details  were  given,  (Ref  4)  for  the 
selection. 

Recall  that  MH  was  defined  as  the  number  of  polar  harmonic 
coefficients,  and  MC  was  the  number  of  control  points.  Some  basic 
guidelines  were  given  for  the  selection  of  MH  (Ref  4:15-18).  The 
value  of  MH  should  be  such  that  MH  <  0.9MC.  The  calculated  coeffi¬ 
cients  should  converge.  The  value  of  the  sidewash  velocity  diverges 
at  a  certain  value  of  MH,  and  the  best  value  of  MH  occurs  just  prior 
to  this  divergence.  The  least  error  occurs  when  MH«MC.  The  run 
time  and  cost  increase  with  the  value  of  MH. 

The  sidewash  criteria  for  selection  was  not  used  because  it 
required  another  computer  program  to  calculate  the  values  for  analysis. 
The  convergence  of  the  coefficients  was  of  no  value  because  the  coef¬ 
ficients  converged  at  all  stations  for  all  values  of  MH.  Different 


values  of  MH  gave  different  coefficients  but  no  analysis  of  the  con¬ 
vergence  could  be  made  because  of  the  random  nature  of  the  numbers. 

The  value  of  MH  =  15  was  selected  by  using  the  example  given  in 
reference  4  and  5,  along  with  analysis  of  the  computer  output  for 
different  values  of  MH.  It  was  noted  that  the  two  examples  in  refer¬ 
ence  4,  where  coefficient  convergence  was  used  as  a  criteria,  gave  MH: 

MH  =  14  for  MC  =  32 
MH  =  14  for  MC  =  24 

For  the  case  of  this  report,  a  value  of  MH  =  14  for  MC  =  30  should  be 
close.  The  sidewash  velocity  criteria  was  used  (Ref  4)  to  give  a  value 
of  MH  =  32  for  MC  =  14. 

The  above  results  indicate  that  the  value  of  MH  should  occur  at 
about  0.5  MC  and  at  a  value  of  MH  just  prior  to  some  kind  of  velocity 
divergence.  If  the  velocity  field  diverges,  the  forces  and  moments  on 
the  store  should  also  show  some  kind  of  divergent  trend.  This  was  the 
approach  taken. 

A  store  location  was  chosen  and  a  trajectory  program  run  was  made 
with  values  of  MH  from  4  through  27.  The  results  were  plotted  for  C^ 
and  C^  versus  MH,  see  figures  13  and  1^.  The  values  of  C^  and  Cffl  both 
show  a  divergent  trend  at  approximately  MH  =  17.  Comparison  with  ex¬ 
perimental  data  revealed  that  C^  matched  best  at  MH  =  10,  while  CN 
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matched  best  at  MH  =  17.  A  value  of  MH  =  15  was  selected  to  give  re¬ 
sults  within  3%  of  experimental  for  both  CN  and  Cm>  The  divergence 
beyond  MH  =  17  was  the  key  to  the  selection,  with  values  of  MH  =  15, 

16,  17  giving  good  results. 

This  method  is  a  very  gross  use  of  the  computer  output  to  select 
a  parameter  within  the  program.  It  uses  a  store  which  is  approximately 
10  feet  long  to  analyze  the  contribution  of  several  fuselage  stations 
to  the  flow  field.  A  better  approach  would  be  to  introduce  a  very 
small  store  and  analyze  the  forces  and  moments  on  it.  This  approach 
would  take  a  large  amount  of  computer  time  and  money,  and  up  to  30 
runs  for  each  station  of  interest. 

The  gross  approach  was  used  at  two  other  store  positions.  The 
results  indicated  new  values  of  MH  for  these  positions.  This  new 
information  was  introduced  into  the  input;  however,  the  final  result, 
with  the  new  values  of  MH,  was  no  more  accurate  than  for  the  run  with 
all  MH  =  15.  The  original  position  analyzed  was  a  critical  position 
in  that  the  solution  at  that  point  gave  the  best  results  for  all  points. 
At  present,  the  selection  of  the  proper  position  for  analysis  can  be 
based  only  on  experience  and  trial  and  error. 

Number  of  XB  Stations  for  Polar  Harmonics 

An  attempt  was  made  to  increase  the  accuracy  of  the  program  in  the 
region  behind  the  wing.  The  number  of  stations  where  polar  harmonic 
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solutions  were  created  was  increased  in  this  region.  Six  extra  stations 
were  added:  however,  no  increase  in  accuracy  was  achieved. 

Free-Stream  Position  of  the  Store 

Confidence  in  the  program  model  of  the  store  was  achieved  by 
moving  the  store  to  a  free-stream  position  10,000  feet  ahead  of  the 
aircraft.  The  values  for  and  Cfn  for  angles  of  attack  of  0°,  2°,  4° 
were  all  within  10%  of  experimental  values. 

Trajectory  Integration 

The  trajectory  option  was  exercised  although  no  wind  tunnel  data 
were  available  for  comparison.  The  stores  were  loaded  on  the  fuselage 
centerline  in  tandem,  with  the  aft  store  separated.  Both  stores  were 
in  positions  which  matched  experimental  positions.  The  program  would 
not  run  until  the  initial  position  of  the  store  being  separated  was 
moved  0.25  feet  below  the  fuselage.  The  error  messages  indicated 
arguments  too  large  for  calculation.  The  program  still  over-predicted 
the  forces  and  moments  when  the  store  was  close  to  the  fuselage. 

The  required  theory  to  predict  store  trajectories  was  established 
in  Chapter  II.  The  required  input  for  the  computer  model  was  generated 
and  the  program  was  used  to  predict  the  forces  and  moments  on  a  separat¬ 
ed  store  at  various  positions  under  the  aircraft  model.  The  results  of 
this  prediction  will  now  be  compared  with  experimental  data  from  wind 
tunnel  tests  on  the  same  configuration. 


Number  of  Polar  Harmonic  Coefficients 


VI.  COMPARISON  WITH  MIND  TUNNEL  RESULTS 


The  comparison  must  begin  with  consideration  of  the  basic  differ¬ 
ences  between  the  actual  model  and  the  Nielsen  Program  model.  The 
program  does  not  model  the  canard  or  tail.  Only  the  wing-fuselage 
combination  is  considered.  The  program  does  not  allow  some  of  the 
fuselage  detail  to  be  modeled.  An  example  is  the  gap  between  the  in¬ 
lets  and  the  fuselage.  The  polar  harmonic  solutions  for  the  near  field 
use  the  exact  fuselage  cross-section,  but  the  number  of  available  sta¬ 
tions  is  limited. 

The  store  model  would  also  be  a  problem  if  the  store  was  other 
than  axi symmetric.  Only  the  aircraft  fuselage  can  be  modeled  with  a 
non-circular  cross-section.  The  store  is  modeled  as  an  equivalent  body 
of  revolution  regardless  of  its  shape.  The  procedure  for  modeling  a 
non-axi symmetric  store  would  be  the  same  procedure  as  described  in 
Chapter  IV  for  the  fuselage. 

Reference  Dimensions 

The  data  comparison  cannot  be  made  until  all  force  and  moment 
coefficients  are  non-dimensional i zed  with  respect  to  the  same  reference 
areas,  or  lengths. 

The  basic  forms  for  the  coefficients  are 


k'  . 


n 


r# 


The  terms  and  SR  were  the  same  for  experimental  and  program  data. 
The  value  of  £R  was  not  the  same.  The  experimental  data  were  based  on 
a  £r  equal  to  the  store  length,  127.5  inches.  For  the  program  £R  is 
equal  to  the  maximum  store  diameter,  15  inches.  The  experimental  data 
were  corrected  for  the  difference  in  £R. 


c  _  Pitching  Moment  127.5  in\ 
Lm  "  q^S^  15.0  in j 


C  = 
m 


original 

Cm 


(83) 


(8.5) 


The  data  for  comparison  were  generated  by  moving  the  store  under 
the  fuselage  with  the  store  CG  as  a  reference.  Table  1  gives  the  store 
position  data. 

Table  1 

Store  Positions 


Run  Sequence  1  Figures  15  and  16  were  the  primary  sequence  for 
comparison.  The  Mach  number  was  0.6  and  the  store  was  positioned  under 
the  fuselage  centerline  with  a  separation  of  approximately  three  store 
diameters.  The  polar  harmonics  of  the  fuselage  were  optimized  for 
this  position. 

The  program  accuracy  in  the  region  of  the  optimum  polar  harmonics 
was  within  ten  percent  of  wind  tunnel  data  with  less  than  3%  error  in 
CN  and  Cm  occurring  at  x=-33  feet.  Cm  showed  less  accuracy  in  the 
area  ahead  of  the  wing.  The  results  diverge  from  experimental  values 
at  x=-51  feet.  This  will  be  discussed  in  a  later  paragraph.  The 
computer  execution  time  for  run  sequence  1  was  27  seconds  on  the  CYBER 
74  system.  This  is  a  typical  value  for  all  run  sequences. 

Run  Sequence  2  Figures  17  through  20  considered  a  y  variation  with 
the  same  z  position  as  run  sequence  1.  The  Mach  number  was  again  0.6. 
The  store  was  positioned  under  the  fuselage,  but  offset  three  feet 
from  centerline. 

The  CN  and  Cm  predictions  were  within  15  percent  of  wind  tunnel 
data  in  the  region  of  polar  harmonics.  The  values  of  Cy  were  accurate, 
but  the  values  were  very  small.  The  yawing  moment,  Cn,  was  underpre¬ 
dicted,  but  this  is  consistent  with  the  near  zero  values  of  Cy. 

The  store  position  is  in  the  region  where  interference  is  the 
greatest,  the  edge  of  the  fuselage.  The  canard,  which  is  not  modeled 


would  affect  the  experimental  data  ahead  of  the  wing,  and  the  wing 

would  create  interference  at  the  wing-fuselage  junction. 

Run  Sequence  3  Figures  21  and  22  were  used  to  consider  a  store 

position  close  to  the  fuselage,  approximately  one  store  diameter 

clearance.  Wind  tunnel  results  were  available  aft  of  x  =  -37.22  feet. 

The  CN  results  were  within  22  percent  of  wind  tunnel  results  until 

x  =  -51  feet,  where  a  spike  in  the  data  occurs.  The  C  results  were 

m 

less  than  50  percent  accurate. 

The  program  accuracy  decreased  when  the  store  was  placed  close 
to  the  fuselage.  This  was  expected  because  of  the  potential  flow  model 
used  by  the  program.  The  model  does  not  consider  any  viscous  effects 
in  a  region  where  the  interference  effects  are  important. 

Run  Sequence  4  Figures  23  and  24  were  made  to  test  for  the  upper 
Mach  number  limit  of  the  program.  The  store  position  was  the  same  as 
in  run  sequence  1.  The  Mach  number  was  0.9,  which  is  above  the  critical 
Mach  number,  but  was  the  only  high  Mach  number  where  data  were  available 
for  comparison. 

A  Mach  number  of  0.9  required  a  new  run  from  the  source  program, 
although  the  same  inlet  velocity  ratio  was  used.  This  induced  an  error 
in  the  results  because  the  flow  through  the  inlet  would  not  be  the  same. 
The  magnitude  of  this  error  could  not  be  estimated,  and  a  complete  new 
model  was  not  generated. 


The  polar  harmonics  were  not  changed.  The  area  of  the  fuselage 
in  the  wing  region  was  optimized  for  a  Mach  number  of  0.6.  This  region 
is  the  region  of  interest;  however,  the  area  upstream  might  not  be 
correctly  represented  due  to  the  unknown  flow  pattern  created  by  the 
higher  than  critical  Mach  number. 

The  predicted  values  for  show  the  trends  of  the  wind  tunnel 
results;  however,  the  accuracy  was  not  as  good  as  the  lower  Mach  number 
case.  The  Cm  results  show  the  proper  trends,  but  the  curve  is  shifted. 
This  shift  occurs  in  the  region  ahead  of  the  wing,  and  is  probably  due 
to  the  unknown  flow  field  caused  by  the  high  Mach  number.  The  program 
was  not  designed  to  handle  this  Mach  number,  but  it  does  give  results 
that  can  be  used  for  comparison. 

The  data  spike  previously  mentioned  occurs  in  all  run  sequences  at 
x  =  -51  feet.  The  wind  tunnel  results  show  an  upward  trend.  This  trend 
is  overpredicted  by  the  program.  The  trend  in  the  experimental  data  is 
a  natural  result  of  some  disturbance  in  the  flow  field  around  the  fuse¬ 
lage.  This  flow  field  is  affected  by  the  components  of  the  aircraft 
and  loaded  store.  The  store  position  is  in  the  downwash  region  of  the 
wing  and  the  region  directly  under  the  horizontal  tail.  Both  are 
possibilities  for  the  cause  of  the  experimental  results;  however,  the 
tail  is  not  modeled  by  the  program.  The  fuselage  is  relatively  wide 
which  may  also  account  for  the  interference.  With  the  flow  around  the 
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fuselage  as  the  possible  cause,  its  effects  would  be  increased  in  the 
program  results  due  to  the  lack  of  a  horizontal  tail. 

An  unsuccessful  attempt  was  made  to  isolate  the  cause  of  the 
spike.  The  program  does  not  model  the  tail;  therefore,  these  effects 
could  not  be  eliminated.  The  effects  of  the  store  which  remained  on 
the  aircraft  were  eliminated  by  removing  the  store  from  the  model. 

The  effects  of  the  inlet  region  on  the  downstream  flow  were  eliminated 
by  using  a  z  direction  traverse  at  various  x  positions  under  the 
fuselage.  Any  disturbance  from  the  inlet  region  should  be  felt  along 
a  line,  at  an  angle,  downstream  of  the  inlet.  This  line  should  be 
predicted  by  the  two  traverses.  No  correlation  could  be  made  from 
the  results. 
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Figure  15.  Run  Sequence  I,  CN,  Mach=0.6,  zRef=6.07  feet 


Fuselage  Station  -  x„  (inches) 


^  ** 


Figure  17.  Run  Sequence  2,  C.,,  Mach=0.6,  yR  *=3  feet 


*•« 


Figure  18.  Run  Sequence  2,  C.  Mach=0.6,  yD  ,=3  feet 

IT1  K6T 


Fuselage  Station  -  xg  (inches) 

Figure  20.  Run  Sequence  2,  C  Mach=0.6,  y_  .=3  feet 


Figure  21.  Run  Sequence  3,  C.,,  Mach=0.6,  zD  ,=3.82  feet 


Fuselage  Station  -  xg  (inches) 

Figure  22.  Run  Sequence  3,  C  ,  Mach=0.6,  zD  -=3.82  feet 


Fuselage  Station  -  xR  (inches) 


Figure  24.  Run  Sequence  4,  C  ,  Mach=0.9,  z„_f=6.07  feet 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  advanced  fighter  design  was  successfully  modeled  using  the 
Nielsen  store  separation  and  trajectory  program.  This  model  was  used 
to  determine  the  loads  on  a  single  store  at  various  positions  under  the 
aircraft  fuselage.  Results  of  the  program  prediction  were  compared 
with  wind  tunnel  values. 

The  computer  model  was  a  potential  flow  model,  with  a  variation 
to  account  for  a  non-circular  crosssection  fuselage.  Only  the  wing, 
fuselage,  and  two  stores  were  modeled;  however,  the  program  could  be 
used  to  analyze  up  to  ten  stores  with  appropriate  pylons  and  racks. 

As  a  result  of  generating  the  computer  model  and  using  it  to  pre¬ 
dict  the  forces  and  moments  on  a  store,  the  following  conclusions  were 
reached: 

1.  The  Nielsen  program  can  predict  the  forces  and  moments 
on  a  separated  store.  For  the  region  of  interest,  where 
the  program  is  optimized,  the  prediction  accuracy  was 
within  ten  percent  of  wind  tunnel  results.  The  region 
of  interest,  at  Mach  number  0.6,  was  under  that  portion 
of  the  fuselage  centerline  which  is  occupied  by  the 
wing.  For  the  same  region  of  interest,  but  Mach  num¬ 
ber  0.9,  the  prediction  accuracy  decreased.  For  this 
case,  results  are  only  within  28%  of  experimental 


2.  The  generation  of  the  program  input  requires  considerable 
knowledge  and  experience.  Flow  regions  of  interest  must 
be  selected  for  optimization.  This  requires  a  knowledge 
of  possible  store  carriage  requirements  as  well  as  store 
free  fall  characteristics.  Parameter  values  within  the 
program  input  must  be  selected  with  very  little  guidance 
given  for  how  the  selection  should  be  made. 

3.  The  program  is  easy  to  use  once  the  model  is  generated. 
Various  store  positions  and  trajectories  can  be  studied 
with  few  changes  required  in  the  input. 

The  Nielsen  approach  to  the  store  trajectory  problem  has  been 
exercised,  and  the  following  areas  were  identified  for  further  study 

1.  The  selection  of  the  optimum  number  of  polar  harmonic 
coefficients  to  describe  the  non-circular  fuselage  is  a 
major  area  of  interest.  The  methods  for  selection  have 
been  discussed  and  an  extension  of  these  methods  was 
used  to  select  the  numbers  for  this  report.  Other 
models  should  be  generated  and  analyzed  to  increase 
confidence  in  the  process. 

2.  In  order  to  build  confidence  in  the  program  it  should 
be  used  to  model  other  aircraft;  however,  caution 
should  be  exercised  when  using  the  program  to  generate 


I 


data  on  an  unknown  aircraft  design.  The  model  should 
be  generated  for  a  low  subcritical  Mach  number.  This 
model  should  then  be  applied  to  a  store  of  known  char¬ 
acteristics  and  at  a  store  position  away  from  the  air¬ 
craft.  The  store  should  then  be  studied  at  numerous 
positions  to  determine  the  limits  of  the  program  at 
that  Mach  number.  Other  Mach  numbers  should  be  ana¬ 
lyzed  in  the  same  manner  to  determine  the  Mach  number 
limits.  This  procedure  should  be  used  to  build  con¬ 
fidence  in  the  program  prior  to  analysis  of  a  store  of 
interest.  For  this  report,  confidence  in  the  program 
was  achieved  by  comparison  with  experimental  results. 

3.  The  effort  required  to  generate  the  model  in'put  could 
be  cut  by  writing  additional  programs.  The  automated 
calculation  of  crossflow  coordinates  and  slopes  for 
the  polar  harmonic  coefficients  would  save  time.  A 
program  to  plot  the  source  generated  shape  and  the  real 
shape  would  help  to  visualize  the  accuracy  of  the  source 
program.  The  programs  would  save  time  and  effort,  but 
not  necessarily  increase  program  accuracy. 

The  Nielsen  approach  to  the  store  separation  and  trajectory  problem 
can  provide  the  information  required  to  determine  what  happens  to  a 
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weapon  when  it  is  released  from  an  aircraft.  The  program  input,  once 
generated,  is  easy  to  use  and  minimum  changes  are  required  to  study 
various  store  positions  and  trajectories. 
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25.  Aircraft  Configuration 


Fuselage  Station 
(inches) 

O 

Total  Area  (in  ) 
(uncorrected) 

Radius  EBR  (in2) 
(corrected) 

522 

5052.1 

36.26 

556 

5131.7 

36.61 

592 

5119.4 

36.56 

632 

4098.8 

31.8 

657 

3821.0 

30.38 

677 

1060.0 

6.84 

Table  3 

Equivalent  Body  of  Revolution  -  Polynomial  Coefficients 


Station  0-75 
75  -  155 
155  -  246 
246  -  327 
327  -  592 
592  -  657 
657  -  682.8 


Cl  C5 ^6 

0.0001497  +  0.355x  -  1.279x2 
0.0098898  +  0.127118X 
0.03729  +  0.007463X 

0.1122185  +  0.662765X  -  0.666058x2 
0.051765  +  0.002558X 
0.1355  -  0.09326X 

1.1766  -  1.1666X 


Ai rcraf t 


Test  Conditions 
Mach  =  0.6 

Total  Temperature  =  100. 3° F 
Total  Pressure  =  1598.2  psi 

=  314.6  lbf/ft2 

Fuselage  angle  of  attack  =  2.5° 
Flight  path  angle  =  0.0® 

tech  =  0.9 


Total  Temperature  =  99.6®F 
Total  Pressure  =  1271.9  psi 

q  =427.1  lbf/ft2 


Run  Sequence  1 


Mach  *  0.6 
yRef  *  0 

=  6.07  feet  below  fuselage  reference  line  (FRL) 


■Ref  (feet) 

CN 

Sn 

Part  Number 

1.51 

0.726 

-0.160 

2393 

9 

0.795 

-0.179 

2387 

15 

0.705 

-0.138 

2381 

21 

0.498 

-0.105 

2375 

27 

0.743 

-0.194 

2369 

33 

0.79 

-0.162 

2363 

39 

0.78 

-0.176 

2354 

45 

0.911 

-0.218 

2348 

51 

1.086 

-0.238 

2342 

54 

0.957 

-0.174 

2339 

Note: 

Part  Number 

refers  to  an 

experimental  run 

number. 
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Run  Sequence  2 


-.W 


Mach  =  0.6 

yRef  =  3-°  feet 

zRe^  =6.07  feet  below  FRL 


,f  (feet) 

h 

c«, 

cy 

V 

c* 

Part  Number 

1.51 

0.771 

-0.165 

0.131 

-0.227 

0.023 

2393 

9 

0.809 

-0.182 

0.13 

-0.201 

0.025 

2389 

15 

0.746 

-0.15 

0.153 

-0.289 

0.024 

2383 

21 

0.550 

-0.112 

0.203 

-0.346 

0.025 

2377 

27 

U.709 

-0.181 

0.095 

-0.034 

0.025 

2371 

33 

0.775 

-0.169 

0.082 

-0.173 

0.022 

2365 

39 

0.784 

-0.177 

0.119 

-0.182 

0.025 

2356 

45 

0.906 

-0.214 

0.066 

-0.04 

0.024 

2350 

51 

1.045 

-0.229 

0.021 

-0.001 

0.023 

2344 

54 

0.924 

-0.172 

0.044 

-0.097 

0.022 

2341 

Run  Sequence  3 

Mach  =0.6 
yRef  =  0 

zRef  =3.82  feet  below  FRL 


xRef  (feet) 


Part  Number 


Run  Sequence  4 

Mach  =0.9 
yRef  =  °*° 


zRef  =  feet  below  FRL 


lf  (feet) 

CN 

C 

m 

Part  Number 

1.51 

.786 

-0.184 

2616 

9 

.898 

-0.228 

2610 

15 

.860 

-0.188 

2604 

21 

.481 

-0.084 

2598 

27 

.751 

-0.216 

2592 

33 

.896 

-0.224 

2586 

39 

.891 

-0.213 

2580 

45 

1.024 

-0.28 

2574 

51 

1.305 

-0.314 

2568 

54 

1.088 

-0.198 

2565 

95 


Store 


=  54 

XX  3,4 

slug-ft 

yy  '  157'5 

slug-ft 

zz  ■  157 • 5 

slug-ft 

xy  ”  *xz 

'yz  '  0 

CG  location  (fus.  Sta.  ft)  =  5.3125  ft 

Length  =  10.625  ft 
Body  Diameter  =  0.625  ft 
Fin  Semi -Span  =  1.45833  ft 


The  Input  deck  for  the  aircraft  will  be  developed.  The  input  for 
the  ogive  store  is  included  for  reference  (figure  27).  The  program 
input  format  is  also  included  (see  figure  28). 

Item  1  NCARDS=3 

Item  2  See  figure  27.  The  data  from  these  cards  is  not 
used  directly  by  the  program. 

Item  3  NSECT=7  1<NSECT<7 

Item  4  NSECT  values  of  the  end  points  of  the  polynomials 
used. 

Item  5  Polynomial  coefficients,  see  figure  27. 

Item  6  (listed  below) 

XSFST  =  .001,  recommended  value,  x/i 

XSLAST  =  1.00357  End  of  body,  including  wake,  x/«. 

XRMAX  =  0.771  point  of  maximum  radius,  x/i. 

XINIT  =  0.01  x/J i 

XFINAL=  1.0  x/ l 

DELX  *  0.02,  recommended  value 

RMAX  =  0.54  maximum  radius,  r/i 

Item  7  NRAT=5,  1<NRAT<5 

This  item  should  not  be  confused  with  item  3.  NRAT  refers 
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to  the  number  of  segments  into  which  the  body  will  be 
divided  for  the  specification  of  the  source  distribution. 
The  segments  need  not  coincide  with  the  NSECT  segments 
except  at  the  end  of  the  body;  the  segments  of  NRAT  must 
cover  the  entire  length  of  the  body  to  be  modeled,  figure 
12. 

Item  3  See  figure  27.  NRAT  values  of  segment  end  points 
Item  9  NRAT  values  of  PERCR.  PERCR  is  input  as  a  fraction 
of  the  local  body  radius  of  the  segment.  A  value 
of  0.9  would  mean  that  the  source  spacing  in  that 
segment  would  be  0.9  times  the  local  radius. 


100 


in  «r  m 

ten 

uuu 

*-»-*- 


sc  r»  ou  9 

E  X  t  t 

UJUbU 


m 

mi 

in 

v£ 

K.  • 

itl 

M 

fh  p4 

• 

O 

•  1 

1 

• 

K>  a  %£  \C  N. 

cm  sc  in  cm  so  in 

K  #  in  K)  OL" 

♦  inhhNN(T>  H  o 

MPNOODO  I  O 

#«HoiJO  IH  I  •  « 

•••••*♦  Ib*  B*  • 


m 

H 

B1 

in  *■* 

so 

Ki 

an 

U  h- 

l>0 

K 

u  u 

3  -J 

• 

• 

• 

<  2 

M 

X  > 

<  SB 

in 

rr  v 

lO 

O  >  Pr 

IT 

o 

m 

O  Uj 

m 

o 

s0 

K  x 

CM 

• 

K1 

Ql  r 

1  X 

• 

• 

III  B  »  tP 

O  • 

k  u  z  m  cm  in 

3(/)U  ♦fT'trcMvO  s© 

*00««*ffc  k  c  N  N  k  in  \fi  U-» 

m  u  i  NorNnHinsn  cr 

mo  o  »i  Hinn  bo  cm 

^  O  O  O  I  O  H  to  CM  • 

•  ••art*  •  H 


Figure  27.  Source  Program  Input 
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Figure  27.  Continued 


Figure  28.  Source  Program  Input  Format 
(Ref  5:  61) 


Figure  28.  Concluded 
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Figure  29.  Source  Program  Output 
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Figure  29.  Continued 
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Trajectory  Program  Input,  Output 


TRAJECTORY  PROGRAM  INPUT 


The  program  input  consists  of  46  items,  figure  30.  Each  item  will 
be  listed;  however,  not  all  definitions  will  be  included.  Reference  5, 
pages  23-40,  must  be  used  along  with  this  report.  The  program  input 
format  is  given  in  figure  31. 

Use  of  the  proper  coordinate  system  is  critical  in  this  input,  and 
numerous  errors  are  possible.  The  basic  reference  frame  is  the  fuselage 
system,  figure  32. 

Item  1  NCARDS=3 

Item  2  See  figure  30.  The  information  from  these  cards 
is  not  used  directly  by  the  program. 

Item  3 

ALFAC=2.5  (degrees),  GAMF=0  (degrees), 

FMACH=0.6,  RH0=0. 0001397  slugs/ft3 

VINF=671.2  feet/sec 

These  values  were  calculated  or  given  in  the  experimental 

data,  see  Appendix  1. 

Item  4  NFU=2,  NPY=0,  NRACK=0 
NSTRS=2  (0<NSTRS<10) 

Item  5 

FLTHC=56.4  feet 


FRMAX=3.05  feet 


Item  6  Lists  the  number  of  sources  used  to  describe  the  EBR. 
This  number  is  printed  at  the  top  of  page  3  of  the 
source  distribution  output,  figure  29,  Appendix  B. 
Item  7  Lists  the  x/A  locations  of  the  sources  used  to  des¬ 
cribe  the  EBR. 

Item  8  Lists  the  non-dimensional  source  strengths.  Item  7 
and  8  are  printed  on  page  3  of  the  source  program 
output,  figure  29,  Appendix  B. 

Items  9,10,11  Apply  only  if  NFU=1  (Ref  5:24) 

Item  12 

NXSTAT=10  NXSTAT<20 
MC=30  MC<100 

VD1NF=0.862  inlet  velocity  ratio 
0.0<VDVINF<1.0 

NXSTAT  is  the  number  of  Xg  station  where  polar  harmonic 
solutions  are  to  be  generated. 

MC  is  the  number  of  control  points  to  be  used,  in  the  range 
of  0°  to  180°  to  obtain  the  polar  harmonics.  The  coordinate 
system  is  shown  in  figure  33.  Note:  the  origin  of  this 
system  is  located  at  the  nose  of  the  actual  fuselage,  see 
figure  32. 


3 


'4 

r-r* 


Item  13  is  repeated  for  each  Xg  station,  NXSTAT  values, 
see  figure  30. 

XSTAT  Xg  location  of  crossflow  plane 
REQ  radius  of  EBR  at  Xg  station 
DRDXC  dr/dx  of  EBR  at  Xg  station 

MPH  number  of  polar  harmonics  to  be  used  at  each  Xg  station 
MPH<50 

NOPSUR  indicates  if  a  portion  of  the  contour  is  open  to  the 
flow,  see  figure  34. 

THOPB  value  of  e  at  which  open  surface  begins 
THOPE  value  of  0  at  which  open  surface  begins 
The  values  of  REQ  and  DRDXC  come  from  the  source  program 
output  for  the  fuselage. 

The  locations  Xg  were  chosen  by  analysis  of  the  fuselage 
cross-section.  The  values  of  Xg  and  reasoning  for  each 
are  given  in  Table  4.  The  values  of  THOPB  and  THOPE  were 
found  during  calculations  for  item  14.  Selection  of  MPH 
is  discussed  in  Chapter  V.  The  value  of  MC=30  was  chosen 
because  it  gave  a  good  range  for  the  possible  values  of 
MPH,  and  it  gave  a  reasonable  representation  of  the  fuse¬ 
lage  cross-section  without  a  large  amount  of  calculations. 

A  computer  program  could  be  written  to  handle  the  calcula¬ 
tions  required  for  items  13  and  14,  but  for  this  case  all 
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calculations  were  done  by  hand. 


Table  4 

Xg  Stations  for  Polar  Harmonics 


Station  (Xg) 

Inches/Feet  _ Reasoning 


-182/-15.17 

highest  point  on  canopy 

-246/-20.5 

start  of  inlet 

-267/-22.25 

center  of  inlet,  open  to  flow 

-297/-24.75 

end  of  inlet  open  to  flow 

-327/— 27 .25 

rapidly  changing  cross-section 

-357/-29.75 

end  of  rapid  change 

-417/-34.75 

point  over  wing 

-487/-40.45 

point  over  wing 

-522/-43.5 

end  of  wing 

-592/-49.33 

point  aft  of  wing 

-598/-49.83 

final  position  so  that  harmonics  at 

station  592  will  be  calculated 

Item  14  is  repeated  for  each  Xg  station  where  data  are  to 
be  input.  The  number  of  cards  is  equal  to  the 
number  of  control  points  MC. 

180 

A9  =  -jqg-  =  6  degrees 

ei  =  -y-  +  ( i -i ) A©  i=l,2,...,MC 

The  value  of  rg  was  determined  using  the  coordinate  system 
of  figure  33  applied  to  the  actual  fuselage  cross-section 
at  each  Xg  station,  figure  35.  The  values  of  Yg  and  Zg 
for  each  control  point  were  determined  from  the  drawing. 
These  coordinates  were  then  converted  into  an  rg(e).  The 
value  of  drg  (e)/de  (feet/radian)  were  determined  from  a 
plot  of  rg  versus  e,  figure  36. 

Analysis  of  the  values  given  in  the  example  (Ref  5: 
125)  revealed  that  the  slope  at  any  point  was  determined 
by  AR/A6  applied  to  the  points  on  either  side  of  the 
desired  point.  This  worked  for  points  where  the  curve 
was  smooth.  At  points  with  sharp  corners,  aR/ao  was 
applied  to  the  desired  point  and  an  adjacent  point,  see 
figure  36.  The  choice  was  made  to  always  pick  the  point 
after  the  desired  point  as  the  adjacent  point.  Judgement 


was  required  as  to  when  the  two  point  method  was  required. 
Item  15 

XBW0C=-28.88  feet  Xg  location  of  the  wing  leading 
edge 

ZBW0=0.333  feet  Zg  location  of  the  wing  leading  edge 
WIOO  incidence  angle  of  the  wing  root  chord  rela¬ 
tive  to  fuselage  Xg  axis. 

Careful  reference  should  be  made  to  the  algebraic  signs  of 
the  input  values,  figure  32.  Values  were  taken  from  the 
engineering  drawing  of  the  wind  tunnel  model. 

Item  16 

CRW=+13.67  (feet),  SSPAN=+16.42  (feet) 

Values  are  from  engineering  drawing. 

Item  17 

NCW=8,  MSW=5  (MSW<30) 

The  values  of  NCW  and  MSW  were  those  recommended  (Ref  5: 
27-28). 

Item  18  numbers  and  specifies  the  location  of  each  wing 
(left)  vortex  trailing  leg,  figure  30.  The  wing 
trailing  and  leading  edge  sweep  angle  as  well  as 
the  wing  dihedral  angle  are  also  specified. 
Guidelines  are  given  (Ref  5:27-28)  for  the  proper 
positions. 


Item  19  is  associated  with  the  wing  twist  and  camber 
NTAC=0  (no  twist  or  camber) 

NUN 1=1  (omit  if  NTAC=0) 

Item  20  is  omitted  if  NTAC=0 

Item  21  specifies  the  wing  thickness  distribution 

NCWS=12  (NCWS  x  MSW<400) 

NUNIS=1  (similar  thickness  distribution  at  all  stations). 
The  value  of  NCWS=12  was  chosen  because  data  for  this  item 
and  item  22  were  available,  see  figure  30. 

Item  22  THETAL(J)  -  slope  of  the  wing  thickness  distribu¬ 
tion  at  the  centers  of  the  thickness  panels.  For 
NUNIS=1  only  data  for  the  chordwise  row  adjacent 
to  the  root  chord  is  input.  The  first  value  is 
for  the  panel  at  the  leading  edge. 

The  procedure  for  finding  THETAL  is  given,  (Ref  5:30-31). 
The  airfoil  for  this  report  was  NACA,  65A005  and  the  data 
for  a  NACA  65A006  airfoil  (Ref  17)  were  scaled  to  provide 
the  required  NACA  65A005  airfoil  section.  This  airfoil 
was  plotted  to  give  the  required  slope. 

Items  23-27  are  associated  with  the  pylon  if  present. 

(Ref  5:31-33) 

Items  28-31  are  associated  with  the  bomb  rack  if  present. 
(Ref  5:33-34). 
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Item  32  is  a  deck  of  cards,  with  one  card  for  each  store, 
figure  30. 

NUMSTR(J)  store  identification  number  <  99 
NSHAPE(J)  store,  shape  identification  number  <  99 
SLTHC(J)  length  of  store,  feet 
SRMAY(J)  maximum  radius  of  store,  feet 
XSNC(J)  location  of  store  nose  relative  to  the  wing 

chord  leading  edge  immediately  above  the  store, 
feet;  positive  ahead. 

YSN(J)  z  location  of  the  store  nose  measured  from  the 

W 

fuselage  centerline,  feet;  positive  to  the  right. 
ZSN(J)  Z  location  of  store  nose 

relative  to  the  wing  chord  leading  edge  immed¬ 
iately  above  the  store,  feet,  positive  below. 
SIC(J)  store  incidence  angle  measured  relative  to  the 
wing  root  chord,  degrees;  positive  nose  up. 

The  store  position  is  input  as  a  function  of  the  store  nose 
with  respect  to  the  wing  of  the  aircraft.  The  incidence 
angle  of  the  store  is  then  input  with  respect  to  incidence 
angle  of  the  wing.  All  experimental  data  were  given  with 
the  store  center  of  gravity  (CG)  located  with  respect  to 
the  store  nose  and  fuselage  reference  line  (FRL).  The 


store  incidence  was  given  with  respect  to  horizontal.  The 
geometry  for  the  store  location  is  given  in  figures  37  and 
38. 

The  store  was  first  located  at  a  reference  position. 
From  this  position  the  store  could  be  moved  without  the 
need  to  go  back  through  the  geometry. 

The  reference  position  was  chosen  to  match  a  position 
where  experimental  data  were  available. 

Xg  =  -33  feet,  Z  =  6.08  feet  (below  FRL),  yB=0,  a=6° 

Store  dimensions  are  shown  in  Appendix  A.  For  this  position 
XSNC(J)  =  +2.192  feet,  ZSN(J)  =  4.256  feet.  The  geometry 
for  ZSN(J)  is  shown  in  figure  38.  The  aircraft  is  at 
angle  of  attack  2.5°,  and  the  wing  incidence  is  0°.  There¬ 
fore,  the  store  incidence  with  respect  to  the  wing  is  3.5°. 

The  calculations  can  be  verified  by  the  trajectory 
program  output  which  gives  the  store  CG  position  with 
respect  to  the  nose  of  the  aircraft. 

The  directions  for  calculation  of  YSN(J)  are  incor¬ 
rect  (Ref  5:34).  YSN(J)  is  measured  from  the  fuselage 
centerline.  This  error  was  discovered  when  the  program 
was  run  with  YSN(J)?<0. 
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Item  33  is  one  card  which  specifies  NSHPT=1,  the  number 


of  different  shapes  for  which  source  distribution 
are  to  be  input. 

Item  34 

MSHAPE=1  shape  number  of  store  to  be  represented 
by  the  following  source  data;  equal  to  one  of  the 
values  of  MSHAPE(J),  item  32. 

MS0R=59  number  of  sources  from  page  3  of  source 
output  for  the  store 
Items  35,36 

DUMX(J)  x/fi.  of  store  sources 
OUMQ(J)  source  strengths 

The  data  for  items  34,  35,  36  come  from  a  source  program  run 
to  model  the  store.  The  procedure  is  identical  to  that  for 
the  fuselage.  The  complete  derivation  for  the  source  input 
is  given  in  Reference  5,  pages  13-13.  The  input  deck  is 
given  in  figure  27  of  this  report. 

Item  37 

NEJECT=2  identification  number  of  the  store  being 
separated 

NSEG=40  number  of  segments  the  body  is  to  be 

broken  into  for  the  force  calculations: 
NSEG  <  40.  It  was  found  that  the  best 


results  occurred  when  NSEG=40. 

NSEGX0=40  number  of  segments  to  the  flow  separation  loca¬ 
tion.  The  recommended  value  of  40  was  used. 
NGAM=0  trajectory  to  simulate  a  wind  tunnel  captive 
store 

NP0LY=2  number  of  polynomials  representing  the  store 

shape.  The  wake  can  carry  no  forces,  and  is 
not  included. 

NR0LL=1  rolling  moment  to  be  calculated 

NFMP=1  empennage  present 

NDAMP=0  no  damping  to  be  included. 

Item  38  (from  Ref  15) 

SMASS=29.8497  slugs  store  mass 

FIXX=5.4  slug-ft3 

FIYY=157.5  slug-ft3 

FIZZ=157.5  slug-ft3 

FIYZa0 

FIXZ=0 

FIYY-0 

Item  39 

XM0M=-5.3125  store  position  where  moments  are 
to  be  taken. 
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XBAR=0  x  positive  forward  of  store  CG,  relative 
to  moment  center 

YBAR=0  y  positive  right  of  store  CG,  relative 
to  moment  center 

ZBARO  z  positive  below  of  store  CG,  relative 
to  moment  center 

For  this  case  the  moment  center  and  CG  are  the  same  point. 

Item  40 

XEND(J)  x/i,  of  end  points  of  the  polynomials 
specifying  the  separated  store. 

Item  41 

C0EF(J,K)  coefficients  of  the  polynomials 
representing  the  separated  store. 

The  data  for  items  40,  and  41  come  directly  from  the  source 

program  input,  figure  27. 

Item  42 

CA=0.34  store  axial-force  coefficient,  reference 
area  is  store  cross-sectional  area.  CA 
was  given  in  the  experimental  data. 

CDC=1.2  crossflow  drag  coefficient, the  recommended 


IPLNR=0  cruciform  empennage 

MSF=5  number  of  spanwise  control  points  on  each  fine 


MSF  must  be  odd,  5<MSF<11.  The  recommended 
value  was  used. 

Item  44 

XTAIL=-9.57  x  location  on  the  store  axis  where  the 
empennage  forces  act  feet,  negative 
RADAV=0.625  average  store  radius  in  empennage  region, 
feet 

FINSS=1. 45833  tail  fin  semispan,  feet,  positive 
PHIR0L=0  initial  fin  orientation  0<PHIR0L<90 
CLALPH=3.1  lift-curve  slope  of  two  exposed  panels 
joined  together,  per  radian.  The 
reference  area  is  the  store  maximum 
cross-sectional  area. 

RADAV,  FIWSS,  and  PHIROL  were  taken  from  the  experimental  data 
Reference  15.  CLALPH  was  determined  using  store  data  and 
DATCOM  (Ref  18:4.1.3.2-4).  The  reference  area  used  in  DATCOM 
is  the  fin  (wing)  surface  area.  This  must  be  redefined  using 
the  store  cross-sectional  area. 


p. 


i  ~ 

! 


t 


i 

3 


» * 

I  4ti 

•J"  J 

I 


F* 


|>';- 


XTAIL  was  determined  using  DATCOM  (Ref  18:4.1.4.2-3) 

which  calculates  the  aerodynamic-center  location. 

Item  45 

This  card  specifies  initial  velocities  of  the  store 
relative  to  the  parent  aircraft.  For  this  case  all 
were  zero. 

Item  46 

DTIME  =0.5  integration  time,  seconds 
TIMEI  =  0.0  initial  time  in  seconds 
TIMEF  =  0.0  final  time,  seconds.  With  TIMEI, 
equal  to  0.0,  the  forces  and  moments  are  given  for  the 
initial  position  of  the  store. 

Item  47  is  used  when  a  trajectory  is  being  restarted. 
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Figure  30.  Trajectory  Program  Input 
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Figure  30.  Continued 
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15.00000 

•  691 C  0 

.15600 

89.93449 

o.ooc 

21.09000 

.71*00 

.24500 

92.06100 

.000 

27.00000 

.74300 

•34100 

92.34723 

.000 

33.30C0C 

.76600 

•50203 

90.43457 

.002 

39.00000 

•84500 

.65300 

90.05194 

•  00  2 

45.00000 

.93700 

•923C0 

90.43125 

O.OOC 

51.00000 

1.04900 

1.30800 

89.72919 

.001 

57.00000 

1.19600 

1.  : 520  0 

89.85390 

•  OC  3 

63.00000 

1.42500 

2.47200 

92.96157 

•  006 

65.00000 

1.73000 

1.58100 

110.13059 

.006 

75.00000 

1.84700 

.70600 

144.08107 

.040 

61.  00000 

1.6740C 

•10500 

167.79309 

.011 

f 7.00000 

1.86900 

••02400 

177.73570 

.007 

53.00000 

1.86*00 

.09600 

180.05=62 

.007 

99.00000 

1.69  >00 

.29200 

180.21292 

.00°. 

105.03000 

1.93200 

.52600 

179.76998 

1.191 

111.00000 

2.00100 

.79000 

179.70391 

1.216T 

117.00000 

2.09SC0 

1.07400 

179.95909 

1.2371 

123.00000 

2.21=00 

1.46400 

179.58490 

.837= 

129.00000 

2.40400 

1.94700 

179.99598 

.503? 

135.30000 

2.62600 

2.68500 

179.36353 

.2541 

141.00000 

2.94500 

2.81830 

It)  7.301 27 

.0495 

147.00000 

3.21500 

1.7*400 

203.52206 

.  014C 

153.00000 

3.31700 

1.98100 

212.15323 

.0031 

155.00000 

3.63100 

2.62000 

211.16546 

-.0156 

165.00000 

3.89600 

2.45200 

222.82844 

-.0375 

171.00000 

4.14600 

1.45500 

241.66195 

-.0841 

177.00000 

4.20200 

•43800 

261.04923 

-.0613 

POLAR  HARMONIC 
COEFFICIENTS 
K  A(K> 

1  -.1S712E*00 

2  .28371E«00 

3  -» 2560  RE *0  0 

4  . 16789E*0 0 

5  -.92550E-01 

6  .5630  4E-0  1 

7  -.41266E-D1 

9  .289558-01 

9  16257E-0  1 

10  . 6869 1E-0 2 

11  -.21291E-02 

12  .47055E-03 

13  -• 70  359E-0  4 

14  .639348-05 

15  -.26714C-06 
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Figure  39.  Continued 
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cm  NO.  3 

IS  POLAR  HARMONIC 

C0ErFICIENTS 

EQUIVALENT  BODY 

theta  range  WHERE 

THE 

XF»  FT 

RADIUS. 

FT  OR/DX 

body  contour  is  nct 

SOLID 

'22.250  00 

2.66000 

• 137B9 

99.300  DEG.  TO  104. 

000  OEG. 

SHAPE  OF  CONTOUR 

THETA 

RADIUS 

OR/DTHETA 

SLOPE  OF  CONTOUR 

0E6. 

FT. 

FT/RAD 

3E6. 

DNU/DX 

3.00000 

.66700 

.00700 

92.35372 

.000  72 

9.00000 

.67400 

•11200 

89.56525 

.000  c  3 

15.00000 

.69100 

•18700 

89.85720 

O.OOOOC 

21.00000 

.71500 

■27000 

90.31234 

o.oococ 

27.00000 

.74500 

.38800 

89.46926 

.00227 

33.0  C  000 

.79500 

.51700 

99.96347 

.00241 

39.00000 

.85400 

.70000 

89.65948 

.00272 

45.0  0000 

.9370C 

.94800 

89.S6565 

.00405 

51.00000 

1.05400 

1.31600 

89.69165 

. 00  26  C 

57.0C000 

1.21200 

1.91300 

89.35671 

.00338 

63.00000 

1.45230 

2.53500 

92.80326 

•003R2 

69.00000 

1.75200 

2.26430 

106.73455 

-.01097 

75.00000 

1.94100 

.71200 

144.65599 

-.04386 

r 1 .0 0000 

1.90C00 

-.26700 

179.99918 

-.013=4 

ET.OCOOO 

1.86600 

-.01500 

179.97389 

.43750 

93.00000 

1.88600 

-.09800 

135.97452 

.82125 

59.00000 

1.90700 

24.79100 

103.39870 

.00900 

105.0  C  COO 

4.63400 

2.24200 

169.18159 

•  036s  7 

111.00000 

4.85600 

2.43700 

174.35006 

.02774 

117.00000 

5.13600 

-7.06900 

260.95962 

.02651 

123.00000 

4.41500 

-6.35300 

263.20274 

.11651 

125.00000 

3.82100 

-4.87300 

273.89942 

.1019' 

135.0  0000 

3.42100 

-2.55400 

265.91029 

.1104* 

141.00000 

3.09900 

1.51000 

205.02216 

.1174* 

147.00000 

3.25C00 

'  1.12400 

217.92225 

.00427 

153.00000 

3.32500 

1.56700 

217.76652 

.00207 

159.00COO 

3.58700 

2.23500 

217.07365 

-.0325' 

165.00000 

3.80100 

1.77600 

229.95598 

-.065*1 

171.00000 

3.95100 

1.0420C 

245.22575 

-.08s 11 

177.00000 

4.02300 

•30400 

262. 67'  63 

-.1014- 

POLAR  HARMONIC 
COEFFICIENTS 
K  A(K) 

1  -.R3303E-02 

2  • 77467E-0 I 

3  -.55936E-01 

4  .  34U7E-01 

5  -.26176C-01 

6  .17893E-01 

7  -.7484SE-02 

8  .30629E-03 

9  . 1 04  35E-0  2 

0  -.74006E-03 

11  .26263C-03 

,2  -.58377E-04 

3  a  B2S15E-Q5 

4  -.68471E-06 

3  a25579E-07 
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Figure  39.  Continued 
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STATION  *C.  4 


XF«  FT 
-24.75000 


THETA 
DEG. 
3.0C000 
9.0  0003 
15.00000 

21.30000 
2  T.O OCOO 
33.30030 

39.30000 
45.0 0000 
51.00000 
57.00000 
63.0 0000 
69.00000 
75.00000 
€1.00000 
€7.00000 
93.00000 
59.0  COOO 

1C5.9C000 
111.00000 
117. C  COOO 
1 2  3.0  C  000 
129.00000 
133.0  9000 
1*1.00000 
147.00000 
153.03000 
159.00000 
165.00000 
171.00000 
177.00000 


15  POLAR  HARHCNIC  C3ESFICIENT$ 


EQUIVALENT  BOOT 
RADIUS.  FT  DR/OX 
2.77000  .07830 


THETA  RANGE  WHEe E  the 
BOOT  CONTOUR  IS  NCT  SCLIO 
0.300  DEG.  TO  C.000  DEG. 


RADIUS 

FT. 

•66-00 
•67630 
.69100 
.71500 
.75300 
.80400 
.86500 
.955C0 
1.06700 
1.23000 
1.47600 
1.69600 
1.795C0 
1.856C0 
3.25400 
4.45600 
4 *565 CO 
4.76293 
4.95300 
5.27700 
5.05300 
4.33700 
3.67700 
3.50900 
3.26400 
3. 332 CO 
3.46700 
3.57400 
3.66300 
3.70500 


r-.' 
U  i 


POLAR  HARMONIC 
COEFFICIENTS 


K 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


A(K) 
-.16694E-01 
•52443E-01 
-.42120E-01 
•  24854E-0  1 
-•  1 2288E-0 1 
. 51 739E-0  2 
-•  18115C-0  2 
•5179 3E-0  3 
-.12269E-03 
•25423E-04 
-•48470E-05 
.61637E-06 
-• 1 0577E-0  6 
•€752RE-08 
-.33577E-09 


SHAPE  OF 
DR/DTHETA 
FT/PAD 
.02300 
•1G9CQ 
.19300 
.29800 
•39700 
•  S380  0 
.74000 
•96500 
1.37500 
1.98500 
2.26200 
1.54900 
.76600 
13.06500 
12.50000 
1.01900 
1.44300 
1.90200 
3.11500 
>2.11300 
-4.33200 
>5.57300 
•3.92400 
-2.90500 
•64800 
.97100 
1.16300 
•93800 
•62100 
.23100 


CONTOUR 

SL0*E  OF  CONTOUR 
OEG. 

91.02902 
89.94032 
89. 39467 
89.06056 
89.20071 
89.21137 
89.45330 
89.11405 
88.81146 
88.76432 
96.12525 
116.59375 
141.89007 
89.39529 
101.59139 
170.11906 
171.45821 
173.22764 
163.B3373 
228.92200 
253.60690 
271.10936 
270.34519 
270.62039 
225.77111 
226.75300 
230.45635 
240.29432 
251.37795 
263.43233 


(e)  Page  5 

Figure  39.  Continued 


ONU/DX 

o.ooooc 

-.00032 
O.OOCCO 
O.OOOOC 
0.00000 
-.00133 
-.00024 
.00022 
.00177 
.00 2- 7 
.00262 
.00440 
.136-6 
.072=1 
.11442 
•09358 
.07353 
•0445s 
•oso3e 

.01337 

.07653 

.07311 

.06-05 

•06°26 

.07-15 

•0C215 

-.03671 

-.07057 

-.08455 

-.05837 
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STATION  NC.  5  15  POUR  HARNONIC  CBErFICIENTS 


EQUIVALENT 

BOOT 

THETA  RANGE  WHERE  THE 

XF«  FT 

RADIUS*  FT 

DR/DX 

80DT  CONTOUR  IS 

NOT  SOLID 

27.25000 

2.97000 

.01930 

0.131  0E6.  TO 

0.000  DEG. 

SHAPE  OF 

CONTOUR 

THETA 

RADIUS 

DR/DTHCT  A 

SL03E  OF  CONTOUR 

DEC. 

FT. 

FT/RAD 

BEG. 

DNU/DX 

3.00000 

.66300 

.02700 

90.68541 

-.00032 

9.0  0000 

.67500 

.10600 

90.07533 

C.OOOCO 

15.00000 

•69100 

.19100 

83.54362 

-8.00608 

21.00000 

.71500 

.23300 

99.06056 

.00089 

27.00000 

.7530C 

.39400 

89.37962 

-.00095 

33.00000 

.79900 

.53600 

83.14477 

-.00106 

39.0  0000 

.86400 

.73700 

86.53551 

o.ocooo 

45.00000 

.95600 

1.01900 

89.17236 

-.00285 

51.00000 

1.07600 

1.42000 

83.15293 

-.00174 

57 .0  COOt 

1.24700 

1.58600 

33.12452 

-.0029= 

t  3.3 0000 

1.49100 

2.31700 

95.76146 

-.0051* 

69.00000 

1.71500 

3.78PQQ 

33.35843 

.020=5 

75.00000 

2.26000 

8.34600 

30.15167 

•02467 

: 1.00000 

3.476C0 

11.72500 

97.50769 

. 06  If  5 

67.00000 

4.72300 

.32700 

173.03941 

•010=5 

93.00000 

4.75600 

•40700 

178. 10577 

0.00030 

99.0  0000 

4.30600 

•75400 

180.06370 

0.00000 

105.00000 

4.91200 

1.42300 

179.84377 

•00092 

111.00000 

5.13400 

1.94300 

180.28894 

-.01766 

117.3 COOC 

5.32200 

1.24500 

193.93331 

-.00436 

123.0  OOOG 

5 .36b  0  0 

-6.S2500 

263.55649 

.00183 

129.900CC 

4.7Cr00 

-5.74800 

269.67424 

.01622 

135.0  0000 

4.18400 

-4.31S00 

270.89306 

.01314 

141.00000 

3.79000 

-3.16000 

270.99633 

.01790 

147.00000 

3.51300 

-2.15800 

268.56198 

.0190* 

153. C 0000 

3. 335  GO 

-.79500 

256.33250 

.01121 

159.00003 

3.34600 

.03400 

248.41761 

-.01696 

165.00003 

3.34600 

.23300 

251.01662 

-.03458 

171.00000 

3.39500 

.23600 

257.02354 

-.07534 

177.0  OOCO 

3.3970  0 

-.01800 

267.30359 

-.08256 

PCLAR  HARMONIC 
COEFFICIENTS 
K  A(K) 

1  -.51833E-01 

2  .32952E-01 

3  -.63105E-02 
A  -.33093E-02 

5  .18753E-0 3 

6  . 3fa 152E-0 2 

7  -• 41 863E-Q  2 

8  .25561E-32 

7  -.10656E-32 

10  .32119E-03 

11  -.70727E-04 

12  .1U75C-Q4 

13  -.12058E-05 

14  .7980SE-07 

15  -.24494E-08 


(f)  Page  6 

Figure  39.  Continued 


ON  NC.  6 

15  R0L4*  MASHCVI 

I  C3£rFI Cl  ENTS 

EQUIVALENT 

BOOT 

THETA  RANGE.  WHERE  THE 

KF,  FT 

RADIUS i  FT 

DR  /OX 

9307  CONTOUR  IS 

NOT  SOLID 

■23.753  00 

3.00  700 

.00256 

0.130  DEG.  TO 

0.000  OEG. 

SHAPE  OF 

contour 

THETA 

RADIUS 

OR/OTHETA 

SLORC  OF  CONTOUR 

OES. 

FT. 

FT/RAO 

OEG. 

ONU/OR 

3.00030 

•66700 

0.00000 

93.30000 

-.00256 

5.0C030 

.67570 

.11000 

89.74427 

-.00268 

15.3GG00 

•  651 C  0 

.16100 

89.54862 

-.0027B 

21.00003 

.714  00 

•  2640  0 

89.41906 

-.00327 

27.0COOO 

.75000 

.39500 

89.22572 

-.00311 

33.0C000 

.79500 

.53300 

89.16051 

-.00346 

39.3  0030 

.66430 

. 70800 

39.66729 

-.00495 

43.0  00  00 

.94300 

.56500 

88.75205 

-.00385 

71.03000 

1.06700 

1.34708 

89.38392 

-.00556 

37.30000 

1.23000 

1.96000 

99.11037 

-.00510 

63.00000 

1.46130 

3.07700 

88.39898 

-.00515 

67.30000 

1.67300 

5.21000 

83.77356 

-.00655 

73.0  3000 

2.35500 

13.52000 

84.34854 

-.00535 

-1.3  0000 

4.15400 

5.74300 

116.87979 

-.03160 

:  7.00000 

4.73  700 

-.00900 

177.10840 

.00016 

9  3.0  0000 

4.75600 

.23600 

190.15923 

-.00016 

79*3 0000 

4.60630 

.75700 

180.34390 

0.00000 

103.30000 

4.91300 

1.30500 

190.13027 

.00077 

111.00000 

5.06000 

1.93600 

190.13904 

.00090 

117.0  0000 

3.30:30 

1.40800 

192.14335 

.00217 

123.33000 

5.37730 

-3.65400 

259.43945 

.00915 

129.00000 

4. 78s  30 

-5.3490  0 

267.16167 

.00534 

13S.JC000 

4.243C0 

-4.36600 

270.83166 

.00390 

191.00000 

3.66300 

-3.18800 

270.53167 

.00420 

147.00000 

3.56300 

-2.36900 

270.47179 

.003*7 

133.0  0000 

3.37500 

-1.40800 

265.64517 

.00251 

159.00000 

3.29300 

-. 72900 

261.48274 

-.007=7 

165.00003 

3.22200 

-.64400 

266.30  309 

-.014=1 

171.00000 

3.15600 

-.40100 

268.23437 

-.01619 

177.00000 

3.136C0 

-.15800 

269.88152 

-.01886 

polar  harnonic 

COEFFICIENTS 
K  A(K> 

1  -.41210E-01 

2  •  1 7G8  8E-0 1 

3  .1T1UE-02 

4  -.68442E-02 

5  .32077E-02 

6  .  72430E-03 

T  -.13679E-C2 
3  .13139E-C2 

7  -.57593E-0  3 

10  .17707E-03 

11  -•  3321  BE  *0  7 

12  .61832E-0S 

13  -.66239C-06 

1*  .43390E-0  7 

13  -.13153C-08 


(g)  Page  7 

Figure  39.  Continued 
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STATION  nc 


.  7  15  POLAR  HARHONIC  COE=FICIEKTS 


KF t  FT 

E0UIVA1 ENT 
RAOIUS.  FT 

BODY 

O8  /OX 

THETA  RANGE  UHE8E 
BOOT  CONTOUR  IS  NCT 

THE 

SOLID 

34.75000 

3.00300 

.00256 

0.900  DEG.  TO  0. 

000  DEG. 

THETA 

DEG. 

RAOIUS 

FT. 

SHAPE  OF 
DR/OTHETA 
FT /RAO 

CONTOUR 

SLD=E  OF  CONTOUR 
OES. 

DNU/OX 

3.00000 

•65100 

.02000 

91.24031 

0.00000 

9.0  0000 

.65%  wO 

.00900 

98.21637 

0.00000 

15.00000 

•67300 

.18400 

89.70887 

O.OOOCO 

21.00000 

•696C0 

.25600 

90.90567 

o.oooco 

27.00000 

.72300 

.33800 

92.09523 

o.oooor 

33.00000 

.76900 

.47500 

9 1.296=8 

o.oooco 

38.00000 

.82430 

.66200 

90.22167 

o.ooooc 

45.00000 

•910C0 

.£  7400 

91.15604 

o.oooor 

Ol.OCOOG 

1.01100 

1.29100 

89.23146 

o.ooooc 

57.00000 

1.17CC0 

1.40400 

96.80557 

C . 0  00  0  0 

£3.0 0000 

1.38600 

3.48600 

84.68225 

o.oocco 

69.00000 

1.75200 

4.60800 

83.81722 

0.00000 

75.00000 

2.34900 

9.92600 

83.31159 

o.oseco 

£1.90000 

3.817C0 

11.58200 

93.24031 

o.oocco 

87.00000 

4.75600 

-.02900 

177.34921 

o.ococo 

53.90000 

4.75500 

•22700 

180.26692 

o.oooco 

99.00000 

4.80600 

.76000 

180.01390 

0.00000 

105.00000 

4 .920  90 

1.32100 

179.97073 

o.oooco 

111.00000 

5.08600 

1.97100 

179.81694 

o.oooco 

117.00000 

5.32200 

1.78500 

189.41=83 

0.00000 

123.09009 

5.46000 

-6.27300 

261.9637= 

0.00000 

124.00000 

4.833C0 

-5.63300 

268.5=604 

o.oooco 

1  35.0  0000 

4. 276  CO 

-4.44300 

271.08390 

o.oocco 

141.00009 

3.89700 

-3.21700 

270.53990 

c.oooro 

147.00000 

3.612C0 

-2.39300 

270.52502 

o.ooooc 

153.00000 

3.39200 

-1.77900 

270.67563 

0.00000 

159.00000 

3.242C0 

-1.27400 

270.45320 

0.00000 

165.00000 

3.12700 

-.86000 

27J. 37753 

o.oooco 

171.00000 

3.05700 

-.49000 

270.10636 

o.ooooc 

177.00000 

3.02100 

-.24500 

271.6364= 

0.00000 

POLAR  HARHONIC 
COEFFICIENTS 
K  A(KI 

1  -.39204E-01 

2  .13867E-01 

3  .45069C-0  2 
A  -. 78955E-02 

5  .33636E-02 

6  .65971E-03 

7  -.17084E-02 

8  • 11554E-0  2 

9  -.4896CC-03 

10  .14501E-33 

11  -.3093SE-04 

12  .46975E-05 

13  -• 48477E-0  6 

14  • 3059 OE-0 7 

15  -.89342E-09 


(h )  Page  8 

Figure  39.  Continued 


STAT ION  NO.  S  IS  POLAR  HARMONIC  CDE"ICIENTS 

EQUIVALENT  BOOT  THETA  RANGE  WHERE  the 


XF»  FT 

RADIUS* 

FT  OR/DX 

BODY  CONTOUR  IS 

NOT  SOLID 

40.58000 

2.99000 

.00256 

0.330  DEG.  *0 

0.000  DEG. 

SHAPE  OF 

CONTOUR 

THETA 

RADIUS 

DR/OTHETA 

SLO’E  OF  CONTOUR 

OEG. 

FT. 

FT/RAD 

DEG. 

ONU/DX 

3.00000 

.65100 

.02000 

91.24031 

-.00027 

9.00000 

.65«00 

.00900 

93.21637 

•00C55 

15.00000 

.67300 

.18400 

89.70887 

-.00079 

21.00000 

.69600 

.25633 

90.80567 

. OC 1 03 

27.00000 

.72900 

.33800 

92.09523 

.00050 

33.00000 

.76900 

.47500 

91.29698 

.00256 

39.00000 

.82400 

.66200 

90.22167 

.00363 

45.00000 

.91000 

.87400 

91.15604 

.00178 

51.00000 

1.01100 

1.28100 

89.28146 

.00560 

57.00000 

1.17000 

1.40400 
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88.31159 
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99.24031 

.03027 
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-.02900 

177.34921 

.00055 
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199.26682 

.01779 
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.76000 

190.01390 

•01732 
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1.32100 

179.97078 
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179.81694 

.01635 
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.00208 

123.00000 
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-6.27 300 
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-.01439 

125.00000 
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-S. 68 300 
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-.01048 

135.00000 

4 . 27i  00 
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-.00532 
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-3.21700 

278.53990 

-.00570 
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-2.39300 

270.52502 

-.00731 

153.00000 
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270.67563 

-.00534 

159. 00000 
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-1.27400 
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-.00916 
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-.00460 
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POLAR  HARMONIC 
COEFFICIENTS 
K  A(K) 

1  -.46237E-01 

2  . 1 761 IE -0  1 

3  .47747C-02 

4  -.S6S46E-02 

5  .44861E-02 

6  .S1305E-03 

7  -.19735C-02 

8  . 1 3977E-0  2 

9  -.60256E-C3 
0  . 18095E-03 

1  -.39015E-04 

2  .59773E-05 

3  -• 62154E-0 6 

4  . 39490E-07 

5  -.U606E-0£ 
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’’C7  STAMON  SIC.  9  15  POLAR  MARNONIC  COEFFICIENTS 


EQUIVALENT 

BODY 

THETA  RANGE  WHERE 

THE 

XF.  FT 

RADIUS*  FT 

DR/DX 

BOOT  CONTOUR  IS  NCT 

SOLID 

43.50000 

3.02200 

.00255 

0.130  0E6.  TO  C. 

000  DEG. 

SHAPE  OF 

CONTOUR 

TMETA 

RADIUS 

dr/dtheta 

SLOPE  OF  CONTOUR 

DEC. 

FT. 

FT /RAD 

DEG. 

DNU/DX 

3.0  3030 

•65000 

.06003 

87.72610 

. OC 273 

9.00000 

•66000 

.09000 

31.23483 

.00272 

15.00000 

.67000 

.20000 

99.37924 

.00263 

21.00000 

•70030 

•30000 

87.90141 

•00126 

27.00000 

.73030 

.38000 

09.50086 

.00243 

33.00000 

•  79C  00 

.51000 

99.82149 

.00230 

39.00000 

.84100 

.63000 

99.04234 

.00203 

45.00000 

•  919  C  3 

.96000 

89.75000 

.00199 

51.00000 

1.04400 

1.33000 

99.13059 

.00051 

57.00000 

1.20700 

1.97000 

66.49540 

.00053 

63.0  0  000 

1.44COO 

3.45000 

85.65523 

.00053 

69.00000 

1.82000 

5.25000 

97.98554 

-•000°5 

75.00000 

2.550  C  0 

11.75000 

fi  7.24453 

-•002c 1 

51.00000 

4.17000 

S.50000 

116.25189 

-.03455 

87.00000 

4.76000 

3.25  000 

142.67579 

-.01360 

53.0  0000 

4.820C0 

•50000 

177.37764 

.02730 

99.03000 

4 .970  C  0 

•77000 

190.01529 

.04337 

105.00000 

4.99000 

1.33000 

180.07573 

.04773 

111.00000 

5.15000 

1.62000 

183.53935 

.03142 

117.00000 

5.33000 

1.15000 

194.82451 

.00537 

123.00000 

5.38030 

-5.45000 

259.37033 

-.01639 

129.0  0000 

4 . 780  CO 

-5.34000 

267.37969 

-.00365 

135.90000 

4.25000 

-4.33000 

278.53421 

141.00000 

3.87030 

-3.05000 

269.24213 

0.00000 

147.00000 

3.56000 

-2.33000 

270.61614 

-.00114 

153.00000 

3.37330 

-1.76000 

270.57605 

-.00122 

159.00000 

3.21000 

-1.24000 

270.12119 

-.00128 

165.00000 

3.11000 
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270.28629 

-.00265 

171.00000 

3.04000 

-.52000 

270.70665 

-.00271 

177.00000 

3.00000 

-.25000 

271.76364 

-.00137 

POLAR  HA="ONIC 
COEFFICIENTS 
K  AIK) 

1  - . 19226E-D 1 

2  .50772E-02 

3  -.37020E-02 

9  .46463E-02 

5  1518  5E-0 2 

6  -• 23542E-0  2 

7  • 34597E-0  2 

8  -.23985E-02 

9  .1081 6E-0 2 

10  -.34286E-03 

11  .77841E-04 

12  -.12489E-04 

13  .13521E-05 

14  -.68948E-07 

15  .26932E-0B 
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ST  AT  I  O^l  NO.  10  15  POLAR  HAORCMIC  COEFFICIENTS 
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86.73*71 
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93.33604 
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92.03624 
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39.0  0000 

.86000 

•67030 

91.97892 
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88.3831? 
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90.71574 
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1 .80030 

4.55000 

90.58397 

0.00000 
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3.73C00 
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.10000 

199.93712 
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5.37000 
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0.00000 

129.00000 
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-3.19  000 

270.49840 

0.00000 

147.00000 

3.57000 

-2.43000 

271.24203 

o.oooco 

153.00000 

3.36000 

-1.76000 

270.64597 

0.00000 

159.00000 

3.20000 

•1.29000 

270.9555? 

0.00000 

185.0  0000 

3.09000 

-.86000 

270.55230 

0.00000 

171.00000 

3.02C30 

-.49000 

270.03107 

0.00000 

177.00000 

2.99000 

-.25000 

271.77950 

o.oooco 
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3  .21067E-01 
A  -.73646E-02 

5  -.22153E-02 

6  • 38650C-0  2 

7  -.22629E-02 

8  . 8224  7E-0  3 

S  -• 20  08SE-3  3 
0  .31725E-04 

,1  -.25102E-C5 

2  -.11319E-06 

3  .S2209E-07 

4  -.53801E-08 

5  .2077 7E -09 
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Figure  39.  Continued 
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,  20.  Abstract  (continued) 

"of  three  store  diameters  away  from  the  aircraft,  is  within  ten  percent  of 
experiment  in  the  region  of  interest,  which  is  the  area  under  that  portion  of 
the  fuselage  occupied  by  the  wing. 
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